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ABSTRACT 

We present a wide area 8x8 kpc), sensitive map of CO (2-1) emission around the nearby starburst 
galaxy M82. Molecular gas extends far beyond the stellar disk, including emission associated with 
the well-known outflow as far as 3 kpc from M82’s midplane. Kinematic signatures of the outflow are 
visible in both the CO and Hi emission: both tracers show a minor axis velocity gradient and together 
they show double peaked profiles, consistent with a hot outflow bounded by a cone made of a mix of 
atomic and molecular gas. Combining our CO and Hi data with observations of the dust continuum, 
we study the changing properties of the cold outflow as it leaves the disk. While H 2 dominates the ISM 
near the disk, the dominant phase of the cool medium changes as it leaves the galaxy and becomes 
mostly atomic after about a kpc. Several arguments suggest that regardless of phase, the mass in 
the cold outflow does not make it far from the disk; the mass flux through surfaces above the disk 
appears to decline with a projected scale length of ~ 1-2 kpc. The cool material must also end up 
distributed over a much wider angle than the hot outflow based on the nearly circular isophotes of 
dust and CO at low intensity and the declining rotation velocities as a function of height from the 
plane. The minor axis of M82 appears so striking at many wavelengths because the interface between 
the hot wind cavity and the cool gas produces Ha, hot dust, PAH emission, and scattered UV light. 

We also show the level at which a face-on version of M82 would be detectable as an outflow based on 
unresolved spectroscopy. Finally, we consider multiple constraints on the CO-to-H 2 conversion factor, 
which must change across the galaxy but appears to be only a factor of ^ 2 lower than the Galactic 
value in the outflow. 


1. INTRODUCTION 

Galactic winds (jChevalier fc Gleg3 119851 : 
iHeckman et al.l ll99Qf) may enrich the circumgalac- 
tic and intergalactic medium and affect the evolution of 
galaxies by moving fuel for future star formation from 
the disk to the inter- or circum-galactic medium. Indeed, 
feedback — often attributed to active galactic nuclei but 
potentially also due to intense star formation— is often 
cited as a necessary element to produce the luminosity 
function observed in the present day galaxy population 
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(e.g.. ICroton et aIll2QQ6f ). Feedback due to star forma¬ 
tion is often specifically invoked as the mechanism by 
which galactic disks self-regulate, making it essential 
to produce commonly observed scaling relation s (e.g., 
[Andrews fc Thompsonll2Qlll : lHopkins et al .1120121 . among 
many others). Enormous abundances of dust and gas, 
much of it enriched, are now observe d outside of galaxies 
via redenning an d absorption (e.g.. [M4nard et alll2QlQl: 
IWerk et all 12013 ) and this reservoir of gas is so large 
(e.g., IWerk et~aL l2Q14l ) that the interplay of the galaxy 
disk with circumgalactic gas and dust must represent 
an essential aspect of galaxy evolution. Even in the case 
where material is driven out of a galaxy disk but does 
not escape the galaxy (a “Galactic fountain” scenario, 
iShapiro fc Eieldl Il976[ ). the redistribution of baryons 
may shape the future evolution of the galaxy. 

As a major mode of feedback from the disk into 
the halo. Galactic winds driven by star formation rep- 
res ent a key asp e ct of this problem (see the review 
by IVeilleux et al.l l2QQ5f ). Indeed, there is observa¬ 
tional evidence that outflows driv en by star forma¬ 
tion are pervasive at both low (e.g.. lArmus et al.l[19^ 
iLehnert fc Heckmai] Il995l: iRnoke^et al.ll2QQ5l ) and high 
redshift (e.g.. iNewman et al junn" Based on newer ob¬ 
servations, there are hints that outflows may often en¬ 
train large masses of molecular gas and dust from the 
star-forming disk, though note that molecular winds are 
not a common result of simulations. These starburst- 
driven molecular winds have been observed in emission 
in stacks of nearby ULIRGs (jChung et al.ll2Qlil) and re¬ 
solved nnai_thn_^sks_of the local sta rbursts including 
M82 (jWalter et NGC 2146 (jTsa.i et, all IMm 
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IKrenke] et allllril). NGC 325fi (l^amot,o et a] .112?)?)^ 


and NGC 253 (iBoiatto et, al.l 1201 3a 'l. Molecular winds 

apparently driven by active galactic nuclei have also been 
seen in many syste ms Te.g.. ICicone et al.ll2Q14[) . includ- 
i ng; Markarian 231 (iFernglio et al.ll2QlQ[ ) and NGC 1266 
(jAlatalo et al.ll2Qlll ). An observationally-motivated pic¬ 
ture is emerging for these AGN-driven molecular winds; 
here we focus on the less well-understood phenomenon 
of starburst-driven winds by studying the most famous 
example. 

Given the potential importance of this phenomenon, 
we know relatively little about the physics of molecular 
winds in starburst galaxies on large scales. In particular, 
the evolution of gas as it moves away from the starburst, 
including the interplay of different components of the in¬ 
terstellar medium remains poorly explored. We know 
that the gas must change phases: most of the gas in 
the starbursts that drive the wind is in the molecular 
phas e, but almost all of the circumgalactic gas is ion¬ 
ized (jWerk et al.ll2Ql4 and references therein), so that 
the cold, neutral component must either change phases 
or fall back to the disk. However, this phase change has 
not been resolved observationally. The fate of molecular 
gas associated with a wind also remains unclear: is gas 
merely removed from the starburst region and then cy¬ 
cled back to a new part of the disk (“a fountain”) or is 
it expelled to spend significant time in the halo (a true 
“outflow”), perhaps even escaping the galaxy? 

Because of its proximity, the proto type edge-on star- 
burst and galactic wind galaxy M8 2 (jLvnds fc Sandag^ 
119631: 10 ’Connell fc Manganol Il978f ) represents an ideal 
target to study a galactic outflow in several phases of 
gas and dust simultaneously. Indeed, for a long time 
M82 hosted the only known exa mple of a molecula r 
outflow driven by star formation (j Walter et al.l l2QQ2[ ) . 
In this paper, we present a new sensitive, wide area 
map ofCOJ = 2 ^ 1 emission around M82 (ED- 
The map shows that molecular gas extends for kilopar- 
secs in every direction beyond the central disk, includ¬ 
ing CO emission coinci dent with the ti dal Hi stream¬ 
ers connec ted to M81 (lYnn et al.l Il99^. the extended 
dust halo (lEngelbracht et al.l I2QQ6I: iKaneda et al.ll2Ql(l 
iRonssel et al.ll2QlQf ). and the outflow visible in Ha, scat¬ 
tered UV, and X-ray emission. It also includes kine¬ 
matic information that allows us to show broad, multi- 
component line profiles in the outflow, slowing rotation 
with increasing vertical displacement, and a clear veloc¬ 
ity gradient along the outflow axis. Combining all of this 
information, we argue that the cold gas around M82 can 
be best described by a “fountain”-type model. The cold 
material confines and funnels the outflow, with the in¬ 
terface giving rise to some of the most spectacular emis¬ 
sion features seen in M82. However, the cold material 
cast out of the disk does not itself travel more than a 
few kpc from the galaxy, on average, instead forming 
a spheroidal halo with characteristic scale length « I- 
2 kpc. The cold gas around M82 is extended compared 
to the central starburst, which is approximately 300 pc in 
size, and the stellar disk, which has vertical scale length 
^ 0.4 kpc. However, the distribution of cold material 
is confined compared to the extent of the circumgalactic 
medium observed in other galaxies via absorption line or 
reddening measurements, ^ 10-100 kpc. 

Our starting framework for the geometry of the M82 


outflo w is the paper by iShopbell fc Bland-HawthornI 
(IT^ . They consider the hot and ionized gas in de¬ 
tail and conclude that the Ha arises from two cones with 
base a few hundred pc and slightly ajar from one an¬ 
other. The Ha arises mostly from the walls of this cone, 
which presumably indicates an encasing medium of cold 
material. Given the relatively low resolution of our data, 
we emphasize the fate of the outflow once it is launched. 
For studies of th e molecular outflow n ear its base, see the 
study of M82 by I Walter et al.l (j2QQ2[ ) and recent demon¬ 
stration ofdenLse_g^assocmted with the base of the out¬ 
flow by iKeplev et al.l (j2QI4l ) . 

Throughout the pap er we adopt a distance of 3.6 Mpc 
(jFreedman et al.l Il99^ so that I" 17.5 pc and I' « 
I kpc. The inclination of the gala xy will often be re l¬ 
evant to our a n alysis and f ollowing iGoker et al.l (j2QI3[ ). 
IWestmoqnettel (|2QQ7[ ). and iGoetz et al.l (jI99Q[ ). we will 
assume that M82’s disk is inclined at 80° and the outflow 
to be perpendicular to the disk, inclined at 10° relative 
to the plane of the sky. Because this inclination is un¬ 
certain and even small changes can have large impact on 
the derived outflow velocities, we also show results for 
±5°. For a thorough review of M82 tha t is beyond the 
scope of this paper, we refer the reader to IWestmoqnettel 

(l2nnl . 

In the appendices, we discuss related issues not appro¬ 
priate for the main text. The first appendix quantifies 
the prospect of finding systems like M82 from unresolved 
spectroscopy by artificially inclining the galaxy. The sec¬ 
ond appendix investigates the CO-to-H 2 conversion fac¬ 
tor and conditions within the molecular gas in detail. 
The third appendix explores some systematics related to 
our treatment of the dust spectral energy distribution, 
including the role of collisional heating. The final ap¬ 
pendix discusses an algorithmic approach for correcting 
single dish maps for contamination by stray light picked 
up the in the “error beam” of the telescope. 

2. OBSERVATIONS 

We used the He terodyne Receiver Array (HERA, 
iSchnster et al.l[2QQ3 ) on the IRAM 30-m telescope to ob¬ 
serve a wide footprint around M82 in GO (2-1) emission. 
The data were obtained between January and March 
2007 and then red uced using the HERAGLES pipeline 
(jLerov et al.l 120091 ). The reduction included baseline fit¬ 
ting using polynomials up to fifth order at both the in¬ 
dividual spectrum and cube stage, rejection of patholog¬ 
ical spectra based on comparing their rms noise to ex¬ 
pectations from the radiometer equation, and construc¬ 
tion and subtraction of OEE spectra from the ends of 
the scan legs. We gridded the data onto a cube with 
a 4" pixel scale using a 17" gridding kernel; combined 
with the native 10.5" beam size of the IRAM 30-m at 
230 GHz, this yielded a map with effective beam size 20". 
This is c oarser than other HERA maps of nearby galax¬ 
ies (e.g., iLerov et al.ll20I^ because in order to cover a 
wide area with good baseline stability, the observations 
did not critically sample the 30-m beam. The channel 
spacing of the WILMA backend used to obtain the data 
is 2.6 km s“^. In the Appendix, we describe an iterative 
correction scheme to account for the effects of the 30- 
m error beam (Kramer, Pehalver, & Greve 2013, IRAM 
memo). We apply the error beam correction to the cube. 
The integrated correction to the flux is ~ 10% and lo- 
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Fig. 1.— Peak CO J = 2 —>■ 1 emission from our IRAM 30-m 
observations of M82 on a high stretch (saturating at 50 mK) and 
20" resolution (grayscale). The figure shows the full area surveyed, 
chosen to cover the outflow, stellar disk, an d the Hi stream er in¬ 
duced by the ongoing interaction with M81 (lYun et al.lll993lb For 
comparison, we plo t the CO J = 1 —)■ 0 interferometer map from 
I Walter et al.l (120021 . in contour), which initially demonstrated the 
molecular outflow in M82. Gray lines indicate our adopted major 
and minor axes. 


cally the effect seldom exceeds 20-30% of the ffux. We 
use a Hann kernel to smooth the cubes and work with a 
final velocity resolution of 5.2 km s“^. Finally, we apply 
the estimated main beam (0.58) and forward (0.92) effi¬ 
ciencies of the 30-m at 230 GHz to convert from the mea¬ 
sured antenna temperature to main beam temperature. 
At this resolution, the cube has rms noise per 5.2 km s“^ 
channel of 7 mK in units of main beam temperature. For 
reference, our beam size and observing frequency imply 
~ 17.4 Jy per Kelvin. 

To highlight the extent and sensitivity of the maps. 
Figured] shows the peak intensity of this CO J = 2 ^ 1 
emission on a high stretch over the full survey area. 
We pl ot the OVRO interferometer map bv lWalter et al.l 
(|2002[) in white contour for reference. Figure [2] shows 
the CO J = 2 ^ 1 integrated intensity and intensity 
weighted mean velocity maps over the region where we 
detect CO. This map has dramatic ally increased field o f 
view relative to the OVRO map of I Walter et al.l (|2QQ2f ) 
but much coarser resolution. It has greatly improved res¬ 
olution and sensitivity compare d to the FCRAO single¬ 
dish map bv lTavlor et al.l (I2QQ1D . It com plements the CO 
J = 1 ^ 0 map of iSalak et al.l (|2Q13[ ). which has com¬ 
parable resolution but much lower sensitivity than our 
HERA map, with « lOx higher rms noise in T^b units 
for matched velocity channels. Similarly, it improves on 
the extent an d sensitivity o f the CO J = 3 ^ 2 map 
presented by IWilson et al.l (|2QI2[) . We will see below 
that comparing the CO J = 3^2, J = 2-^I, and 
J = I ^ 0 lines yields interesting constraints on physi¬ 
cal conditions. 


3. DATA AT OTHER WAVELENGTHS 

In order to place the molecular superstructure of the 
galaxy in context, we compare our new CO map to ob¬ 
servations of M82 at other wavelengths. As tracers of the 
ionized wind, we use a narrow-band ima ge of Ha eniission 
from the Local Volume Legacy survey (|Dale et ani2QQ9l: 
iLee et ani2QQ9D: Ha velocity in formation from longslit 
sp ectra bvlMcKeith et al.l (II995 ) and a Fabry-Perot cube 
bv iShopbell fc Bland-Hawthorn (|I998F W e also use far- 
ultraviolet (FUV) emission f rom GALEX (|Hoopes et al.l 
I2QQ5I: iGil de Paz et al.ll2QQ7D. and X-ray emission from 
Chan dra ( Strickland et al.ll2QQ4 iStrickland fc Heckman! 
I2QQ7D . T o trace starlight, we u se the 2MASS K band im¬ 
age from iJarrett et al.l (j2QQ3f ) ; the more sensitive IRAC 
3.6/im map from LVL is contaminated by PAH emission 
in the region of the outfiow. 

To study the distribution of dust in the halo of M82, 
we use infrared o bservations from H erschel at 70, 160, 
250, and 350/im (jRpussel et al.ll20Tol ) and Spitzer maps 
at 3. 6, 8, and 70Rm from the Local Volume Legacy snr- 
(Engelbracht et all 120061: [Dale et al.l 120091: iLee et all 


vey 


12009 ). As discussed in iRpussel et al.l (j20I0[ ). the presence 
of the bright starburst in the field means that studying 
M82’s dusty halo in the i nfrared require s high dynamic 
range imaging. Eollowing iRpussel et al.l (j20I0f )the Her¬ 
schel data have been “cleaned” in a manner similar to ra¬ 
dio interferometry data, though in detail the procedure is 
refined from that paper. The algorithm reconstructs the 
image as a series of point sources, in each case replacing 
the complex Herschel point spread function, which has 
substantial extended structure, with a Gaussian. To es¬ 
timate the Herschel PSE, we combined observations of 
point sources, bright sources for the wings and fainter 
sources for the core to avoid saturation. These mod¬ 
els took into account the rotation of the telescope. The 
IR images were then cleaned by subtracting a scaled ver¬ 
sion of the PSE, corrected for orientation of the telescope 
from the brightest point source in the image and replac¬ 
ing it with a Gaussian containing equivalent fiux. We 
cleaned to a limiting brightness of O.I, 0.04, 0.6, and 0.2 
Jy arcsec“^ at 70, 160, 250, and 350/im. The beam of 
the Gaussian used to replace the PSE at each wavelength 
was: 4.2" (70/im), 8.6" (I60/im), 13.5" (250/im), and 19" 
(350/im), in each case corresponding to about 75% of the 
EWHM of the Herschel PSE. The final “cleaned” Her¬ 
schel data have approximately Gaussian beams for their 
bright emission and we convolve these to a common reso¬ 
lution of ~ 24" for most analysis. Eor most of this paper 
the key point of this analysis is we are able to distin¬ 
guish faint emission even in the presence of M82’s bright 
starburst. The final SPIRE images had a small positive 
zero-level offset, which we fit and subtracted away from 
the galaxy. The effect of this subtraction is shown in the 
appendix. 

Using the matched-resolution, cleaned infrared inten¬ 
sity maps, we build a point-by-point spectr al energy dis- 
tribution, which we fit using the models o flDraine fc Li 

e foll owing the parameter ization of iDraine et al. 

and lAniano et al.l (j2QI2D . The resulting fits in¬ 
clude estimates of the dust surface density, and 

the mean radiation field {U) for each locatiorru Other 


The models distinguish a PDR-like component and minimum 
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Fig. 2.— CO J = 2 —>■ 1 emission from M82 at 20" resolution. The top panel shows integrated intensity, with contours beginning at 
0.2 K km s“^ and increasing by factors of \/2. The bottom panel shows the emission-weighted mean velocity, on a linear stretch from 
80 km s“^ to 330 km s“^. To indicate scale, we plot a black bar of size 1 kpc at the distance of M82. Both panels show the bright, rotating 
central disk surrounded by a region of CO emission that extends for many kpc in each direction. Some of this extended emission may be 
associated with an extended disk, some of it appears associated with the tidal streamers that link M82 to M81 in Hi, and some of the 
material is coincident with the starburst-driven outflow visible in Ha and X-ray emission. 
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parameters in the model, for example the PAH mass frac¬ 
tion, gpAH, and the breakdown of the radiation field dis¬ 
tribution, are less certain and we do not consider them 
in th i s paper (seelEngelbracht et al']l2QQ6l: iKaneda et al.l 

I2Q1QI: iBeirao et al.ll2Q15L for disc^ sion of the PAH distri¬ 

bution). The iDraine fc: Li] (|2QQ7l) models consider only 
radiative heating of dust. Shocks and processing of 
small du st grains are known to be important around 
M82 (e.g. IKaneda et al~]l2QlQl:IBeirao et al.ll2Q15f ). so col- 
lisional heating of dust might become important in some 
areas around the galaxy. However, in the appendix we 
show that for reasonable limits to the bulk gas tempera¬ 
ture and even a modest radiation field we do not expect 
collisional heating to be dominant over most of the area 
around M82. However, more detailed modeling of the 
interface between the very hot, rare wind and the cold 
material is certainly n eeded. The append ix also shows 
that the choice of the IDraine et al.l (|2QQ7[ ) models does 
not substantially bias our results; most of our analysis 
requires a linear tracer of dust mass and we sh ow that 
within a factor of ~ 2, the IDraine et al ] (l2nnl models 
match the results of fitting a modified blackbody mod¬ 
ulo a linear scaling. The appendix also compares the {U) 
from the model fit to the best-fit dust temperature. 

In order to help constrain the physical conditions in the 
molecular gas, we use ma ps of other CO tran sitions. We 
take CO J = 1 —> 0 from lWalter et all (j2QQ2f ). who com¬ 
bined OVRO and IRAM 30-m data and maps of the CO 
J = 3 2 line from the JC MT Nearby Galaxies Legacy 

Survey (|Wilson et al.|[2QT^ . The appendix presents a 
detailed comparison of the three maps, while the main 
text summarizes the results. 

We also compare to Hi 21-cm emission, for which we 
us e several data se ts. The D-conf i gnrati on data set 
of lYnn et ahl (j 19931 ) and lYnn et aP (jl994l ) provides a 
low resolution 70"), high surface b rightness s e nsitiv - 
ity view. Our own reduction of the lYnn et ahl (jl993f ) 
C-configuration data provides a higher resolution view 
(which we convolve to 20" resolution to match the CO 
data). Finally the best matched resolution column den¬ 
sity estimate of Hi comes from our combination of the 
archival VLA data, Effelsberg single dish data, and mod¬ 
eling to convert the Hi absorption into emission. This 
best-estimate column density map has resolution « 15" 
and we convolve it to the 20" resolution of the CO data 
for comparison. We derive intensity weighted mean ve¬ 
locity fields and velocity dispersions at 20" and « 70" 
resolution. We treat these velocity measurements as un¬ 
defined near the disk of the galaxy where Hi goes into 
absorption against the continuum. 

Collectively, these data give us information on the ex¬ 
citation of CO, mostly near the disk; the location and 
kinematics of the outflow as seen in ionized gas (Ha); 
the distribution of small dust grains, likely polycyclic 
aromatic hydrocarbons (PAHs, 8/im); the distribution of 
larger dust grains (70-350/im); the interstellar radiation 
field; the distribution of stars, and the location and kine¬ 
matics of atomic gas (21-cm). We mention additional 
processing as we present individual comparisons. 

4. A WIDE, SENSITIVE CO J = 2 ^ 1 MAP 

Our CO J = 2 ^ 1 observations, shown in Figure [21 
offer the most complete view of the molecular medium 
around M82 to date. Together, the integrated inten- 


Near-IR (2MASS K Band at 20”) 
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Fig. 3.— Our CO (2-1) map (red contours) relative to a high- 
stretch plot of the AT-band light from stars (grayscale) with major 
and minor axes shown by gray lines. CO emission extends well 
beyond the stellar distribution along the minor axis of the galaxy. 

sity map and velocity field show a bright, quickly ro¬ 
tating disk embedded in a distribution of fainter CO 
emission that stretches for kiloparsecs. Indeed, it is 
striking that the columnated morphology familiar from 
UV, X-ray, and Ha images of the M82 wind would be 
hard to identify from the integrated CO intensity map 
alone. Instead, CO emission extends in every direc¬ 
tion, simi lar to results from infrped imaging of dust 
emission JE^ elbracht et al. ] I2QQ6I: iRoussel et al.l I2Q1QI: 
IKaneda et al.ll2QlQlb 

The picture becomes clearer when we compare the CO 
to the stellar disk in Figure [3l CO extends far above 
the stellar disk, with the largest extensions above and 
below the central starburst region along the minor axis of 
the galaxy. This is suggestive of M82’s familiar outflow, 
especially as seen in Ha. The molecular emission above 
and below the disk is not confined to the minor axis, 
however, with CO emission extending above and below 
all along the major axis of the galaxy. 

Even though the integrated intensity is surprisingly 
amorphous, the velocity field shows an impressive de¬ 
gree of structure. The central disk, which is by far the 
brightest feature in the integrated intensity image (note 
the logarithmic stretch in Eigure[2j), shows strong, clear 
rotation. Above and below this disk, a velocity gradi¬ 
ent is visible along the minor axis, a point we will will 
return to below {^. Even above and below the main 
disk, rotation is still visible as a velocity gradient along 
the major axis, but the strength of the rotation seems to 
weaken as one moves away from the main disk. 

This more detailed structure is also visible in the in¬ 
dividual channel maps of CO emission (EigureH)), where 
is it possible to pick out clumps of CO emission far 
above and below the main disk, filamentary structures 
that extend up and down in individual channels (e.g., 
see ulsr = 235, 110, and 79 km s“^), and an ex¬ 
tended low surface brightness region around the main 
disk that includes an eastward extension of the galaxy. 
The channel maps also show that the emission can have 
a very broad velocity distribution. Eeatures often persist 
across several channels and the step between channels 
is Av ~ 15 km s“^. Some of the broad line widths re¬ 
flect rotation in an edge-on disk. However, the broadest 
features tend to be concentrated along the minor axis. 
Eigure[5]shows this by plotting the line of sight velocity 
dispersion calculated from the second moment. 
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Fig. 4.— Individual channel maps of CO emission over the same area as Figure [2] The figure shows every third channel over the velocity 
range of significant emission with the numbers in each map indicating the mean LSR velocity for that 5.2 km s“^ wide channel. A gray 
contour shows emission at 21 mK, or S/N = 3 in an individual channel. Dotted lines show the major and minor axis. A bar in each 
panel shows a length corresponding to 1 kpc for reference and a white circle shows the beam. 
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Fig. 5.— Line of sight velocity dispersion calculated from our CO map. The broadest line widths are concentrated in, above, and below 
the starburst near the minor axis of the galaxy. Some of these broad features show evidence for multiple components, signatures of gas 
surrounding the conical outflow. 
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5. ASSOCIATION WITH THE SUPERWIND AND 
DECOMPOSITION INTO REGIONS 

Figure [6] shows the CO distribution in contours rel¬ 
ative to signatures of the well-known hot outflow from 
M82: Ra (a mixture of in situ and scattered emission), 
scattered far-ultraviolet (FUV) emission, and soft X-ray 
emission from hot gas. Above and below the disk, we see 
qualitative correspondence between the CO and the out¬ 
flow tracers, particularly the Ra. The coincidence might 
be expected given that the FUV light shows the location 
of dust grains scattering emission from the starburst and 
some of that dust is likely to be mixed with molecular 
gas. Similarly, Ra emission traces recombinations, which 
will occur in regions of high gas density, often along the 
surface of neutral clouds. 

In the ne xt sections, we show a kinematic association of 
Hi and CO. I Walter et al.l ()2QQ2f ) have shown a similar as¬ 
sociation of CO and Ra. This suggests that the ionized, 
atomic, and molecular phases are mixed along the minor 
axis. We also see that the neutral (atomic and molecular 
combined) gas matches the morphology of dust surface 
density very well. There is also a good correspondence 
of CO emission with the filaments of PAH emission seen 
in the mid-infrared. Thus the three “cold” phases of the 
ISM — H 2 , Hi, and dust — appear closely associated. 
Modeling the dust spectral energy distribution shows ev¬ 
idence of high interstellar radiation fields along the mi¬ 
nor axis — i.e., the dust temperature is high along the 
outflow — providing additional evidence that this cold 
gas is indeed interfacing with the outflow (though shocks 
or collisional heating at the interface with the hot gas 
might also contribute some to the heating). Moreover, 
individual CO and the Hi spectra along the minor axis 
appear broad and sometimes double-peaked, indicative 
of a multi-directional outflow almost aligned with the 
plane of the sky. 

This evidence suggests that the extended CO emission 
along the minor axis is associated with M82’s superwind. 
Our maps also include the bright starburst region and 
show a large molecular disk that extends along the major 
axis of the galaxy out into the tidal streamers create d by 
the galaxy’s interaction with M81 (jYun et al.lll994l ). In 
order to study the molecular component of the outflow 
quantitatively, we have decomposed the area around M82 
into regions dominated by the outflow, the bright central 
disk, and the extended disk. 

The last panel in Figure [6] shows our region defini¬ 
tion. Based on the CO intensity and velocity, we de¬ 
fine a bright disk region (red). We define the outflow 
region to cover areas of bright Ha, FUV, and X-ray 
emission outside this region (blue; shown in all pan¬ 
els). The remaining area corresponds to the extended 
disk and tidal streamers (gray). The extent of the “out¬ 
flow” (blue) region parallel to the major axis is 3 kpc 
(i.e., ±1.5 kpc about the major axis). For simplicity, we 
do not treat the north and s outh differently, although 
IShopbell fc Bland-HawthornI (jl998f ) have shown asym¬ 
metry between the two in both orientation and opening 
angle. Figure [6] and subsequent figures show that our 
simple region definition captures signatures of hot gas, 
hot dust, the minor axis gradient, and wide line profiles. 
We also indicate a cone with a base of 300 pc and an 
opening angle of 20° (±10°); we will find some evidence 


TABLE I 

Integrated CO Emission 


Region 

CO Flux 
(K km s-i 

CO Luminosity 
(K km s“^ pc^) 

Entire HERA Map 

4.3 X 10® 

5.3 X 10® 

... Central Disk 

3.1 X 10® 

3.8 X 10® 

... Outflow Region 

0.8 X 10® 

1.0 X 10® 

... Outflow 2 :| > 1.5 kpc 

0.12 X 10® 

0.15 X 10® 

... Disk & Streamers 

0.4 X 10® 

0.5 X 10® 


Note. — CO J = 2 ^ 1 emission. Luminosities assume d = 
3.6 Mpc. For region definition see Figure [6] 


for this opening angle from line splitting later in the pa- 
per and the geometry corresponds a pproximately to the 
IShopbell fc Bland-HawthornI (jl998f ) picture of the out¬ 
flow. 

A caveat to this decomposition is that some emission 
above and below the midplane in the outflow region may 
still come from an extended disk. While it is highly in¬ 
clined (i ^ 80°), M82 is not perfectly edge on. We do 
observe an extended disk of CO along the major axis. 
At i = 80°, a 5 kpc disk would be foreshortened to have 
extent of about a kpc along the minor axis, which may 
contaminate measurements of the outflow. We will see 
below that material beyond about 1.5 kpc from the mid¬ 
plane (the two dashed lines in Figure [6j) is very securely 
associated with the outflow. Closer to the disk we show 
evidence that the outflow appears important to the CO 
emission, but expect contamination from a foreshortened 
extended disk. 

Table [T] reports the integrated emission for our whole 
map and individual regions. We report both the flux, 
from summing Iqo n and the luminosity, Ico x U x 
(P. We also note the flux above a height of 1.5 kpc in 
the outflow region, which we can securely assign to the 
outflow. 

6. CO AND OTHER PHASES OF THE COOL ISM 

Neutral atomic hydrogen (Hi) and dust both also ex¬ 
tend to great distances around M82 and represent key 
points of comparison to understand the evolution of the 
molecular phase of the superwind. 

6.1. Atomic Gas 

Although most of the gas in the disk of M82 appears 
to be molecular, the extended tidal features are mostly 
atomi c, so that an enormous Hi supers t ructu re surrounds 
M82 (jCottrelll 119771: lYun et ahl 1 1 9931 Il994l) . Figure [3 
compares the column density and kinematics of atomic 
gas to our CO map. In the top left panel, we see that 
the shapes of CO and Hi features near the starburst 
are similar, but that the Hi features extend well be¬ 
yond CO. The high column Hi arcs to the northeast 
and northwest are the well-kno wn tidal features induced 
by the interaction with M81 (jYun et al.lll994l) . CO is 
present in these tidal features, but as they extend pro¬ 
gressively further from the disk the amount of molec¬ 
ular gas diminishes. The top right panel in Figure [7] 
shows that directly by plotting /h 2 , the fraction of gas 
mass that is molecular (as opposed to atomic); here we 
use aco = 2 Mq pc“^ (K km s“^)“^ for simplicity. Gas 
near the starburst and in the stellar disk is mostly molec¬ 
ular but becomes progressively less so with increasing 
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Fig. 6.— Our CO (2-1) map (red contours) plotted over tracers of the M82 outflow in grayscale: (top left) Ha, (top right) far ultraviolet, 
and (bottom left) soft X-ray emission. The bottom right panel shows the two-dimensional region assignments that we use to break the 
map into emission mainly from M82’s disk (pink), emission associated with the outflow (blue), and emission not clearly associated with 
either — mostly the extended disk and molecular gas associated with the tidal streamers (gray). The outflow region is also shown in as 
a rectangular blue outline in the other three panels. We use this region, which has width dzl.5 kpc about the minor axis, in a number of 
subsequent calculations. The flnal panel shows the major and minor axes and includes lines ±1.5 kpc above and below the major axis. 
The flnal panel also shows a cone with opening angle 20°, which is the opening angle implied by the line splitting that we observe. 


distance from the burst and disk. 

Figure 0 shows that the Hi exhibits the same minor 
axis velocity gradient as the CO. The Hi minor a xis ve - 
locity gradient was first recognized bv iCottrehl (jl977[ ). 
who placed a lower limit of 23 km s“^ kpc“^ on the gra¬ 
dient and noted the link to the minor axis Ha filaments. 
In general, the CO and Hi velocity fields agree. Compar¬ 
ing CO and Hi velocities across the survey area, we find 
no mean offset, though there is a substantial 20 km s“^ 
rms scatter between individual local Hi and CO mean 
velocities. We defer a more detailed investigation to fu¬ 
ture work focused on improved Hi data. Thus from Fig¬ 
ure 0 we conclude that the Hi and H 2 appear mixed 
where both are present, showing similar morphology and 
kinematics, but that Hi is much more extended than H 2 
traced by CO. 

Finally, the bottom right panel in Figure 0 shows that 
broad Hi lines, indicated by a high rms line of sight ve¬ 
locity dispersion, tend to be associated with the minor 
axis of the galaxy, a signature of the complex line profiles 
associated with the outflow. 

6.2. Dust 

We expect dust to be mixed with both molecu¬ 
lar and atomic gas. We also expect that most of 
the dust associated with the M82 outflow will have 
been entrained along with the cooler material because 
the h ot portion of the superwind will q uickly destroy 
dust. iShopbell fc Bland-HawthornI (j 19981 ) estimate the 


timescale for grain destruction by sputtering in the very 
hot (X-ray) portion of the wind to be only ^ 10^ yr (see 
lOstriker fc Silkl[l973l:lDraine fc Salpeterl[l97^ . Dust as- 
sociated with the cooler phase can survive much longer. 
IShopbell fc Bland-Hawthoirnl (jl998f ) also point out that 
dust should even survive for a long time (^10^ yr) in the 
lower temperature ionized gas associated with Ha emis¬ 
sion or a galactic halo, T ^ 10^ K (instead of the 10^ K 
associated with the hot X-ray emitting gas). Therefore 
although dust directly associated with the hot wind will 
be destroyed, if material ejected in a cool phase becomes 
ionized as it leaves the galaxy we do not expect its dust 
to be immediately destroyed. Dust therefore offers the 
potential to uniformly measure gas mass in all but the 
hottest phase of the outflow. 

Figure [8] shows the results of our IR SED fitting, the 
distribution of dust mass surface density, Edust, and 
the mean interstellar radiation field heating the dust, 
{U) (equivalently dust temperature, which might include 
some contribution from collisional heating; see the ap¬ 
pendix). The dust mass distribution shows an excellent 
morphological match to the CO, though low surface den¬ 
sity dust extends further from the galaxy than CO emis¬ 
sion. These extensions all match features in the Hi maps 
(Figure [7]). Together, the three maps show the neutral 
gas becoming more atomic and less molecular at greater 
distances from M82 and both phases mixed with dust. 
Most of the visible dust emission can be associated with 
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Fig. 7.— Atomic gas over the area of our CO survey, (top left) Hi column density (grayscale) at matched resolution with CO contours 
(red contours). Blue rectangles indicate the area of the outflow (same as Figure 0. {top right) Approximate H 2 mass fraction, /h 2 = 
Mh 2/(^H2 + ^Hl), with contours indicating where the gas is 10, 25, and 50% molecular, (bottom left) Hi velocity held at lower resolution 
(beam in the lower left corn er) to take advantage of increased signal to noise. Again the area of the outflow is indicated. Note the minor 
axis gradient (|Cottrelllll977l h similar to what we observe in the CO (Figure [2]). (bottom right) Hi velocity dispersion, also at low resolution, 
tracing the line width over regions of high Hi column. Regions associated with the outflow appear to have higher velocity dispersions, 
reflecting more complex, multi-component line profiles. 


either Hi or CO, so Figure [8] does not show evidence for 
a substantial component of dust associated with ionized 
or hot gas except perhaps above ~ 3 kpc to the north of 
the galaxy. 

The outflo w is particularly evide nt in (U) (right panel). 
As found bv iRoussel et al.l (|2Q1Q[ ). the infrared colors in 
the wind region indicate hotter dust, and thus a more 
intense interstellar radiation held along the outflow. UV 
emission scattered off dust shows that radiation from the 
starb urst escapes along the cavity blown b y the super¬ 
wind (|Hoopes et al .11200^ : I Coker et al.ll2Q13f ). This same 
radiation heats the dust in or along the edge of the cavity. 
Dust in the extended disk is shielded by large amounts 
of intervening material. Inasmuch as dust appears well- 
mixed with gas, this provides further evidence that the 
cool gas along the minor axis is associated with the out¬ 
flow. 

The dusty outflow from the galaxy is particularly 
visible in the 8/im emission, which is d ominated by 
PAH emission (see the spectra l mapping oflBeirao et 51 
IMl)- The Sfim emission (jEngelbracht et al.l l2QQ6^ 


sh ows a filamentary s tructure similar to the Ra (e.g., 
see iKaneda et al.ll 2 QTQh . extending several kpc above and 
below the galaxy. We show the comparison in Figure 
[H where we plot CO contours from individual channel 
maps over the 8 /im emission. CO features in individual 
channel maps show good correspondence to the 8 /im fila¬ 
ments, which in turn show detailed structure down to the 
resolution limit of IRAC, ^ 2 " 35 pc at the distance 

of M82. Until high resolution CO data become available 
over a wide held, the IRAC emission may give the best 
det ailed view of cold IS M structure in the outflow (e.g., 
see IVeilleux et al . 11200 ^ . 

6.3. Molecular Mass Estimates and Dust-to-Gas Ratio 

We wish to translate CO intensity into molecular mass, 
which requ ires an estimate of t he CO-to-H 2 conver- 
sion factor. iBolatto et al.l (|2Q13b[ ) review this topic and 
iBolatto et al.l (|2Q13af ) have considered the CO-to-H 2 con¬ 
version factor for the case of the molecular wind in NGC 
253. In the appendix, we constrain aco in M82 by 
comparing low-J CO transitions and combining Sdust, 
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Fig. 8.— Dust mass surface density (top) and the mean interstel lar radiation field i ll uminating the dust (bottom), derived from point- 
by-point fitting of the infrared spectral energy distribution following IDraine &: Lil (|2007l b [Praine et al.l (120071 b Before fitting, all data were 
convolved to the limiting resolution of Herschel at 350p,m. In the left panel the red contours show CO integrated intensity contours. In 
both panels, blue outlines show our “outflow” region (FigureO. The dust mass contours begin at log^^Q S^ust [M© pc“^] = —I and step by 

0.15. The {U) contours begin at I (a Solar Neighborhood field) and step by \/2. The dashed region shows where we carry out the SED fit, 
which corresponds to an area above a limiting surface brightness of /250 = 15 MJy sr“^. The dust mass distribution shows a remarkable 
match to the CO emission. The faint extensions of Sdust beyond the CO are real and correspond to features seen in the Hi map, so that 
dust around M82 appears mixed with a combination of atomic and molecular gas (see the appendix for an analysis of the dust-to-gas ratio 
and conversion factor based on these data). The right panel shows that the radiation field (equivalently dust temperature) returned by our 
SED fitting shows hotter dust in the region of the outflow. This argues for a direct association between the cool gas and the superwind and 
agree s with the observation of FUV and Ha light scattered off dust grains along the minor axis (e.g., see lHoopes et al.]l2005l : (Coker et al.l 
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Fig. 9. — Contours for individual c hannel maps of CO emission in black contour over the Sfim emission (mostly from PAHs, 
lEngelbracht et ^120061 : IDale et "311 120091') in the the same area as Figure O A blue outline shows our two dimensional “outflow” re¬ 
gion. Individual CO features exhibit a good match to the filamentary structures seen in the mid-IR dust map at much higher resolution. 
The black contours show where peak CO emission across a seven channel range exceeds 21 mK, or S/N = 3 in an individual channel. 
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CO, and H i emis sion. Following IWeifi et ^ (|2QQ1[) and 
IWeifi et all (j2QQ5f ) we conclude that the observed line ra¬ 
tios do not require optically thin gas and indeed appear 
difficult to reproduce for optically thin CO emission. The 
three lowest J ^^CO transitions could arise from either 
hotter, low density (T > 30 K, n < 10^ cm“^) gas or 
cooler dense gas (T < 20 K, n ^ 10^ cm~^ ). Previous 
result s usi ng also using COisotopologues by IWeifi et ^ 
(j200l[) and IWeifi et al.l (j2005f ) argue for the former case. 
This hot, low density gas implies aco ^ 2-3 times lower 
than Galactic, but with significant uncertainties. This 
ag rees with an a lysis o f the virial masses of resolved clouds 
bv iKeto et al.l (j2005f ). 

We use dust emission as an independent tracer of the 
total cool gas to derive another constraint on aco- In 
the appendix, we show that bright CO emitting regions 
in M82 show a high ratio of CO (2-1) to 350/im contin¬ 
uum emission when compared to a large reference sample 
compiled from the HERACLES and KINGEISH surveys. 
This argues for aco ~ 2.5 times lower in M82 than in a 
typical disk galaxy. This is more than an order of mag¬ 
nitude higher than the optically thin value. 

We also use the modeled du st surface density , Ednst 
and follow a modified v ersion of iSandstrom et ffil (j2Q13l ) 
and iLerov et~al] (j2Qllf ) to compare dust, CO, and Hi. 
This approach describes regions of a galaxy using a fixed 
gas-to-dust ratio, and a fixed CO-to-H 2 conversion 

factor, (aco 5 so that 

^^GDR X Edust = <acO X Jco + ^HI • (1) 

Here Edust and Ehi are the mass surface densities of dust 
and atomic gas inferred from IR and 21-cm observations. 
Jco is the CO intensity. Observations of a wide range of 
CO, Hi, and dust values simultaneously constrain the 
gas-to-dust ratio and CO-to-H 2 conversi on factor (see 
iLerov et al.l I2Q11I: ISandstrom et al.l I2Q13L for more de¬ 
tails). Adopting this approach, we simultaneously solve 
for the best-fit aco and 6gbr in the extended disk of 
M82. We find the same result for aco implied by the 
simple /co 2 -i /^350 treatment and a dust to gas mass 
ratio of (^^dr ^ 0-009. 

Our M82 dust -to-gas ratio is 35% lower than the 
value found by ISandstrom et al.l (j2Q13f ). They stud¬ 
ied the CO-bright parts of many disk galaxies using 
the same CO line and dust models and found a dust- 
to-gas ratio of 0.014. Our va lue does resemble the 
0.01 found bv iDraine et al.l (j2007[ ) considering integrated 
galaxy SEDs using the same dust models. Thus M82 ap¬ 
pears somewhat dust-poor compared to the inner part 
of a large disk galaxy. This could agre e with its small 
mass and the mass-metallicity relation (jTremonti et al.l 
120041) and indeed the re is some evidence of snbso- 
lar m etallicity in M82 (jOriglia et al.l l2004 iNagao et al.l 
I 2 OIID despite its apparently n ormal abundance patterns 
(jEorster Schreiber et al.l[200If ). Note that we emphasize 
comparative stat ements when discu ssing the dust-to-gas 
ratio because the lDraine fc Lil (j2007[ ) dust models appear 
to have a normalization that is ~ 2 times too high. Thus 
our best-estimate true dust-to-gas mass ratio would be 
^O.OOb l-to-200). 

Synthesizing (aco estimates from Edust, IR intensity, 
and CO lines, our best estimate (aco is given in Equa¬ 
tion IB8I It is (aco = 2.5 Mq pc“^ (K km s“^)“^in 


the extended disk and faint parts of the outflow 
and (aco = 1 Mq pc“^ (K km s“^)“^in the bright 
part of the outflow and the starburst. We also re¬ 
port results for the simpler case of a single a^^ ^ 
2 Mq pc“^ (K km s“^)“^. Given that our estimates 
span (a^Q^ ^ 1-4 Mq pc“^ (K km s“^)“^, a factor of 
^ 2 represents a reasonable, realistic systematic uncer¬ 
tainty on the conversion to molecular mass. 

These conversion factors, temperatures, and densities 
differ from the Milky Way but are not truly extreme. 
Based on this analysis, one might speculate that the ex¬ 
tended CO emission comes from a close analog of the re¬ 
solved cloud population seen bv iKeto et al.l (j2QQ5[ ) in the 
starburst region of M82. Our knowledge of CO substruc¬ 
ture is circumstantial, but the small size of the bright 
PAH filaments, the “normal” CO line ratios, and the 
similar /co 2 -i /^350 ratios all argue that physical con¬ 
ditions inside the gas may not change overwhelmingly 
between the inner disk and the brighter p a rts of the out¬ 
flow. Reinforcing this idea, iKeplev et al.l (j2Q14l ) showed 
some evidence of dense gas traced by HCO+ extending off 
the starburst disk into the outflow. This would support 
the idea that the material in the M82 outflow consists of 
dense, opaque clouds of molecular gas. 

7. MINOR AXIS MASS PROFILE: CHANGING PHASE AND 
FALLING PROFILES 

With H 2 , Hi, and dust estimates we can examine the 
composition and total amount of cold interstellar mate¬ 
rial along the outflow of M82. EigurefTOl shows our best 
estimate of this profile along the outflow region. The y- 
axis shows mass per unit vertical length integrated over 
a width ±1.5 kpc around the major axis. The x-axis 
shows vertical displacement from the disk with positive 
indicating north. We show profiles for H 2 , Hi, and to¬ 
tal gas estimated from dust assuming a dust to gas ratio 
^^DGR = 0.009. The vertical error bars indicate the sys¬ 
tematic uncertainty in the mass profile. Eor H 2 this is 
dominated by the CO-to-H 2 conversion factor. Eor dust 
this is set by our adopted dust to gas ratio. 

Eigure [TT] shows the corresponding plot for Ha emis¬ 
sion and near-infrared light, which tracks the distribu¬ 
tion of stellar mass. Here we do not estimate mass, but 
plot intensity averaged as a function of vertical displace¬ 
ment. Again we average within ±1.5 kpc of the minor 
axis (the blue “outflow”) region. We normalize the pro¬ 
files to unity near ±0.5 kpc to show the relative decline 
of GO, starlight, and Ha. 

Eigure [To] shows that Hi and H2 both contribute a large 
fraction of the gas mass around M82 between 0.5 and 
1.5 kpc from the plane. H 2 likely makes up most of the 
gas within ^ 1 kpc, though this depends on the adopted 
(aco- In this same regime, the gas implied by dust is a 
good match to the total Hl±H 2 gas. This is formally 
by construction (see appendix) but our res ults would be 
l argely the same if we simply adopted the iDraine et al.l 
(j2007( ) dust-to-gas ratio. 

At larger heights, the surface density of all phases 
declines with increasing distance from the disk. The 
strongest decline is in the molecular component, which 
declines sharply with increasi^ distance to both the 
north and south of the galax30 The associated expo- 

This is not a sensitivity effect; we detect GO over most lines 
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Fig. 10. — Linear mass density profile along the direction of the outflow of molecular gas (blue), atomic gas (red), and total gas estimated 
from dust (purple). The y-axis reports the mass per unit length integrated within ±1.5 kpc of the minor axis (the blue rectangle in Figure 
[3 as a function of height above the plane. That is, the curves are proportional to average projected mass surface surface density and the 
value can be multiplied by length along the minor axis to recover a mass. Here “gas from dust” indicates dust surface density multiplied by 
the estimated ratio of gas mass to dust mass, S^ust x ^dgr with (5dgr = 0.009 (see appendix). Error bars indicate systematic uncertainties, 
e.g., from the CO-to-H 2 conversion factor, dust-to-gas ratio, or calibration. Although we show vertical displacements as low as 500 pc, 
we caution that close to the disk, the association of material with outflow is less straightforward. 


TABLE 2 

Scale Lengths Along the Minor 
(Outflow) Axis 


Quantity 

Scale Length (kpc) 

molecular gas 

0.5 ±0.1 

Ho: intensity 

0.6 ±0.1 

atomic gas 

1.7 ±0.2 

/h2 = H2/(HI + Ha) 

0.9 

dust mass 

1.4 ±0.3 

starlight 

0.4 ±0.1 


Note. — Over the range 0.5-4 kpc from the 
disk. See Figure ITOl The ± values indicate 
the spread between north and south. For an 
outflow velocity of 450 km s~^, it takes about 
2.2 Myr to travel 1 kpc. 


nential scale length is ~ 0.5 kpc and almost the same to 
the north (« 0.42 kpc) and south 0.53 kpc) of the 
galaxy. 

The fall-off in CO closely matches the decline in Ha 
intensities (Figure [TT]). Thus, Figure [6] and the match in 
the profiles in Figure [TT] support the idea that bright 
Ha emission arises from the interface of the cavity 
carved by the hot wind with cold, denser material (e.g . , 
iShopbell fc Bland-Hawthornll19981: IWestmoquettell2007 ). 
As the dense, confining material falls off, the Ha inten¬ 
sity also falls off almost in lock-step. 

Figure [TT] also shows the decline in the stellar profile at 
large heights. The slight vertical misalignment between 

of sight out to at least ±2 kpc, by which time the decline is clearly 
evident. 
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Fig. 11. — Intensity of iF-band light (starlight), Ho; emission, and CO along the minor (outflow) axis. All three profiles are normalized 
to have an average value of unity at ±0.5 kpc. Note the slight misalignment between the CO midplane and the stellar disk midplane 
(also visible in Figure [3]). Ho: and CO are clearly more extended than the stellar disk and decline at a very similar rate at greater offsets 
from the disk midplane. 


the bright CO disk and the stellar disk appears in the 
asymmetry of the north and south profiles (this is also 
visible in Figure [3]). With a scale length of ~ 0.4 kpc, 
the stars have the sharpest decline of any component, 
emphasizing the degree to which the gas — even the 
H 2 — really has been ejected from the central region. 
In particular, CO and Ha emission outside ^ 1.5 kpc 
is very clearly extraplanar and we shall see that this is 
the emission that can best be linked directly to the hot 
outflow. 

The other two neutral gas tracers decline more slowly 
than the CO. For jzj < 3.5 kpc, the Hi shows a scale 
length ^1.7 kpc, that is roughly three times longer than 
the H 2 . This longer scale length of Hi compared to H 2 
means that the phase of the neutral ISM is changing 
systematically as one moves away from the disk. The 
molecular fraction /h 2 = H 2 /(HI ± H 2 ) declines almost 
exponentially with a scale length 0.9 kpc as a func¬ 
tion of vertical displacement from the disk. As a re¬ 
sult, molecular gas contributes substantially to the ISM 
around M82 within about 1.5 kpc of the disk but is sub¬ 
dominant compared to Hi by a vertical displacement of 
3 kpc. For the most part, these statements are robust to 
the conversion factor prescription adopted, though near 
the disk this matters to the exact accounting of where 
H 2 exceeds Hi. 

Dust shows a scale length intermediate between the CO 
and Hi, « 1.4 kpc — steeper to the south and somewhat 


shallower to the north. The dust profile combined with 
fcoR ~ 0.009 is a reasonable match to the total gas 
inferred from adding the Hi and H 2 estimates (the gray 
line). There is not clear evidence for a large amount of 
ionized gas associated with dust around M82. If present, 
this would appear as dust without associated gas traced 
by either CO or Hi. The best candidate for such emission 
would be above 3 kpc to the north of the galaxy, but this 
is near the edge of our maps and faint, so we hesitate to 
ascribe much significance to it. 

We have argued based on dust temperatures. Hi veloc¬ 
ity dispersions, and the apparent kinematic association 
of CO and Hi that the Hi in Figure [10] is probably linked 
to the disk and outflow rather than superimposed high 
latitude material. In this case, given an outflow velocity 
a displacement along the minor axis translates linearly to 
the time since gas was expelled from the starburst and 
the X axis in the profiles can be thought of as a time¬ 
line. If we make the assumption that the starburst has 
continuously expelled gas at a constant rate then we can 
consider this profile as a plot of the time evolution of 
gas as it leaves the starburst. Thus our first-order in¬ 
terpretation of the vertical profiles is that the phase of 
the gas changes as it moves away from the galaxy and 
that gas and dust disappear from the minor axis outflow 
as they reach larger and larger distances from the disk. 
The relevant timescales are all of order a Myr; for an 
outflow velocity of 450 km s~^ (Section l9.1[) it takes gas 
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TABLE 3 

Scale Lengths Over a Wide Field 


Band 

Minor Axis 
[kpc] 

Major Axis 
[kpc] 

70/i,m 

0.8 

0.7 

I60/i,m 

1.1 

0.9 

250/i,m 

1.7 

1.5 

350)am 

1.6 

1.4 

CO 

0.5 

0.6 

Hi 

3.0 

2.6 

AT-band 

0.4 

0.7 


Note. — Fits to data outside 1 kpc. 

“Minor axis” refers to an opening angle 
of 90° about the north and south out¬ 
flow axis. “Major axis” refers to the rest 
of the area. Figure EU shows the rele¬ 
vant geometry and the Figure fT^ shows 
the profiles used for the fit. 

^2.3 Myr to travel a kpc. Therefore roughly every Myr, 
there is an e-folding in the amount of H 2 associated with 
the outflow and roughly every 2 Myr the gas becomes a 
factor of e more atomic. Meanwhile the overall amount 
of gas in the outflow diminishes with an e-folding time 
of ^ 3 Myr. 

For the phase change, the logical explanation is that 
photodissociation by light from the starburst and disso¬ 
ciation in shocks outpace H 2 formation as clouds move 
away from the galaxy. It seems fair to assume that the 
density and shielding of molecular gas will decline com¬ 
pared to the starburst if a cloud is placed in the lower 
ambient pressure of the outflow with less nearby ma¬ 
terial to shield it and no larger molecular medium to 
replenish the cloud mass. Meanwhile the dust model¬ 
ing and scattered UV light both clearly show that much 
of the cool ISM along the outflow is being bathed in 
a strong radiation field that presumably includes dis¬ 
sociating photons. Shocks also appear to play an im¬ 
portant role in the outflow, at the very l east process¬ 
ing s mall dust grains in the outflow (e.g., iBeirao et al.l 
l2QI5f ) and he ating a small fraction o f the H 2 to very high 
temperatures IVeilleux et al.l (|2QQ9[) . The expected rate 
at which /h 2 will change depends on cloud properties 
and geometry (combined shielding and density structure 
and incident radiation field), but we note that in sim- 
ulatio ns treating the opposit e problem (molecule forma¬ 
tion) , iGlover fc ClarkI (|2QI2f ) see dramatic changes in the 
molecular fraction of initially atomic clouds in roughly a 
Myr. It seems reasonable to consider that if one tipped 
the balance in the other direction, a cloud could be dis¬ 
sociated on a similar timescale. 

8 . THE WIDE-FIELD VIEW OF DUST AND GAS: THE 
FINITE EXTENT OF THE COOL MATERIAL 

To help quantify the ultimate fate of cold material 
that is ejected from M82, we also calculate profiles of 
gas and dust emission over all azimuthal angles out to 
large scales. Following the logic of the previous section, 
if there is less material in each successive shell around the 
galaxy then continuity arguments suggest that material 
must stall or fall back. The minor axis profile suggests 
that this stalling or fallback happens. By considering all 
azimuthal angles in this analysis, we account for the pos¬ 
sibility that emission may move out of the outflow region 
that we have defined. Recall that even if the material 


becomes ionized, we would expect it to remain visible in 
dust emission unless the ionized gas is very hot. 

Because we work with faint emission, we use the IR 
images themselves for these profiles rather than the dust 
SED fits. The SED fitting requires significant detections 
over individual lines of sight in many bands, while the 
images allow us to average data into radial profiles and so 
even recover signal at radii where individual lines of sight 
are not detected. This also means that our results are 
sensitive to the zero point of the map, but Eigure[T2]and 
the profiles do not show evidence for either faint extended 
“negative” emission or a plateau to a background level; 
we show the effect of uncertainty in the zero point in the 
appendix. 

We moved all maps to the resolution and astrometry 
of the cleaned SPIRE 350/im map. We then divide the 
galaxy into quadrants oriented along the major and mi¬ 
nor axis. The dotted lines in Eigure [12] illustrate the 
geometry, showing the division into quadrants with the 
circle indicating a radius of 10 kpc. We constructed ma¬ 
jor and minor axis profiles out to very large radius, which 
we plot in Eigure [131 The figure shows average profiles 
for the major (dotted) and minor (outflow, solid) axis 
for each IR band, CO, Hi, and K band. Eor comparison, 
we plot the profiles expected for spherical, optically thin 
distributions with density profiles p oc and a = — 2, 
—3, and —4. The a = —2 profile offers a key point of 
comparison because it shows the borderline between a 
diverging and a converging mass profile. 

Eigures [T^ and [T3l show that the cold material is con¬ 
centrated around M82. The observed IR profiles decline 
steeply compared to a p oc profile. Even extrapo¬ 
lating the observed distribution to infinity beyond the 
sensitivity limits of the data, most of the material would 
still be concentrated near M82. Only the Hi appears 
to follow a profile that implies equal mass in successive 
shells (p oc r“^). Eigure [T2] shows that this is because 
progressively larger circles encompass more and more of 
the extended tidal material. A secondary effect may be 
that most of the gas near the galaxy itself is molecular 
rather than atomic. 

Thus, high-stretch infrared images of a wide field 
around M82 (Eigure [12]) do not suggest that material 
flows indefinitely out into the halo. Inste ad. Eigure [13] 
ygues that the deep Herschel images by iRoussel et al.l 
(|2Q1Q[ ) already capture most of the cold material present 
around M82 and that this material is mainly confined to 
a spherical region a few kpc in size around the galaxy. 

In detail, Eigures [12] and [13] show a number of other 
points important to interpret the cold outflow. Eirst, 
the dust temperature is hotter along the outflow than 
outside the outflow. It also become s cooler as one move s 
further from the galaxy (as shown in iRoussel et al.ir2QlQ[) . 
both along the outflow and outside the outflow. This 
can be seen in the quicker decline of the 70pm and even 
the 160pm bands relative to the SPIRE bands. However, 
while the dust becomes cooler, it does not become so cold 
that the 250pm-to-350pm color changes in an extreme 
way (see the very similar decline of the gray and purple 
lines in Eigure [T3|) . Taking this to mean simply that 
the peak of the blackbody spectrum (in Fj,) drops below 
A = 160pm but not A = 250pm, then the dust cools but 
not too far below ^ 20 K. Dust is not vanishing from the 
IR images because it becomes too cool to see (see also the 
























16 


Leroy et al. 


50'0" 


40'0" 


69“30'0" 




CO > 0.5 K km s 
HI > 10^° cm"' 


SO'O" 


40'0" - 


Maj'ijr:), 


I SPIRE ZSO/xrrTf^ . / 

4 , -■ _r^. ^ Minor 


o '. 


eg^so'O" - 


09''58"'0® 56"'0® 


54"'0® 


PACS 160/xm 


Fig. 12 .— Infrared images at 160/rm (left) and 250/rm (right) shown at high saturation with CO (red) and Hi (blue) contours showing 
the extent of bright molecular and atomic gas. Dotted black lines show the major and minor axis definitions used for the profiles in Figure 
[la which are also labeled in the right panel. The circle indicates a radius of 10 kpc. 



Fig. 13.— Radial profiles of infrared intensity and gas tracers over a wide area around M82. The a^-axis shows displacement on the sky 
from the galaxy center. The y-axis shows intensity normalized to unity at r = 1 kpc. We build separate profiles along the major axis (dotted 
lines) and along the minor axis (solid lines). Gray lines show the projected surface brightness profiles of several density distributions. The 
profiles show that: (1) CO emission declines more rapidly than any other ISM tracer, indicating a changing phase of the gas around M82; 
and unlike any other tracer it is more extended along the major axis, refiecting M82’s extended gas disk (2) dust emission is slightly more 
elongated along the minor axis, likely as a result of the outflow, but this effect is not large; (3) dust temperature, indicated by infrared 
color, declines steadily moving away from the galaxy; (4) Hi is the most extended circumgalactic cool component, likely because it reffects 
the tidal interaction in addition to the outffow; and (5) aside from Hi, all of the other tracers are concentrated towards the galaxy. That is, 
their profiles decline much faster than a projected p on r~‘^ density profile and so there is less emission (and thus mass) in each successive 
shell. Table [2] reports scale lengths associated with these profiles. 
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comparison to mm observations in iRoussel et al.|[2QlQ[) . 

Second, the d ust distribution is ro und on large scales. 
Again following iRoussel et al.l (|2Q1Q[) , who noted a large 
amount of dust not associated with the superwind, note 
that while the IR scale lengths do tend to be slightly 
shorter along the major axis (into the plane of the 
galaxy), this effect is modest, a difference of usually only 
^0.1 kpc or so. This gives rise to the round appearance 
of the highly saturated IR images in Figure [T21 We ex¬ 
perimented with fitting ellipses to the 160/im and 250/im 
isophotes and found axis ratios of ~ 1.1 (ellipticity ~ 0.1) 
for the low intensity contours that encompass a large part 
of the dust distribution. 

Third, the differences between the extended Hi dis¬ 
tribution and the dust emission imply a varying dust- 
to-gas ratio in the cool gas around M82. Much of the 
extended Hi emission at the relatively high contour of 
A^hi = 10^^ cm“^ 1 Mq pc“^) does not have a 
clear analog in the high-stretch SPIRE maps. Some spa- 
tial filtering rn ay be expected from the data processing 
(|Roussell 12013F but a basic comparison shows that we 
would expect to detect 250/im emission from much of 
the extended Hi, especially after averaging, if the same 
dust were mixed with Hi throughout. Much of the Hi su¬ 
perstructure may thus represent gas from the extended 
Hi disk of M82 that was originally poor in dust. 

Fourth, these figures show the confined nature of the 
CO very clearly. Molecular gas is, as we have seen, ex¬ 
tended compared to the stellar disk along the minor axis. 
However, it is confined compared to the other tracers of 
the cool gas. In contrast to the other gas tracers but 
in agreement with the stars, the molecular gas extends 
more along the major axis than the minor axis. 

Finally, the figures shows a quantitative meaning for 
the fr equently quoted extent of ^ 6 kpc for M82’s dusty 
halo (|Engelbracht et aP I2QQ6I: iRoussd et al.l |2Q1Q[) : by 
this distance the cool dust intensity drops by two orders 
of magnitude relative to the intensity at r ^ 1 kpc. At 
the same time, this distance represents about the point at 
which the intensity of dust emission approaches the sen¬ 
sitivity of Herschel. If the decline is exponential, then 
all but a few percent of the mass outside 1 kpc lies in¬ 
side 6 kpc. If decline is closer to p oc then there 
will be equal mass per logarithmic interval of radius, so 
that we expect the same amount of material between 10 
and 100 kpc as between 1 and 10 kpc. In either case, 
the present-day flow of material through successive sur¬ 
faces around the galaxy must diminish with increasing 
distance. 

A major caveat to these mass continuity arguments, 
particularly at large scale, is the assumption that the 
starburst has a continuous, fixed mass outflow rate. Mass 
continuity must apply, but the outflow rate from the star- 
burst itself may vary with the star formation history plus 
some delay time for the relevant feedback processes to 
play out. A natural next step to understand the M82 
outflow will be an analysis that combines the star forma¬ 
tion history of the burst with the radial distribution of 
outflowing cold material. 

9. KINEMATICS OF THE OUTFLOW: VELOCITY 
GRADIENT AND LINE SPLITTING 

So far we have examined the distribution of different 
mass components around M82. Our observations also 


TABLE 4 

Minor Axis (Outflow) Velocity Profile 


Quantity 

Value 

Mean velocity 

... systemic 

211 km s-^ 

... 2 ; < —1.5 kpc 

118 km s-i 

... 2 ; > 1.5 kpc 

269 km s-i 

Mean outflow velocity 

• • • (T’out) sinz 

75 km s“^ 

• • • (L’out) for sin 2 = 10° 

430 km s“^ 

• • • (L’out) for sin 2 = 5° 

860 km s“^ 

• • • (L’out) for sin 2 = 15° 

290 km s“^ 

95% envelope (relative to Vsys) 

... for 2 ; < —1.5 kpc 

138 km s-i 

... for 2 : > 1.5 kpc 

132 km s“^ 

Minor axis gradient (on sky) 

... within ±1.5 kpc of minor axis 

32 km s ^ kpc ^ 

Typical line width 

... at peak for 2 : < —1.5 kpc 

110 km s ^ 

... at peak for 2 ; > 1.5 kpc 

140 km s“^ 

... integrated for 2 ; < —1.5 kpc 

160 km s“^ 

... integrated for 2 ; > 1.5 kpc 

230 km s“^ 


Note. — See Figures \ 2 \ and fT4] 


yield high spectral resolution measurements of the kine¬ 
matics of the gas. These offer a powerful tool to inves¬ 
tigate the mechanics of the outflow. Most obviously, the 
minor axis of M82 shows a ve locity gradient (no ticed in 
the molecular gas as early as iNakai et al.l 119871 ) , which 
is strikingly evident in our velocity field (Eigure [2j) and 
dat a cube. This gradient, which is also present in the 
Hi (ICnttre]lll197^ . lets us estimate the magnitude of the 
outflow velocity and explore the possibility of gas ac¬ 
celerating away from the plane. The line profile along 
individual lines of sight also offers key information. Its 
shape has often been invoked as a way to explore the 
three dimensional geometry of an outflow. Specifically, 
the appearance of two components along an indivdiual 
line of sight (“line splitting”) has b een taken as evidence 
for a tilted conical geometry (e.g. , iMcKeith et al.lflMsI: 
iShopbell fc Bland-Hawthor3ll998[ ). 

9.1. Position-Velocity Diagram, Minor Axis Gradient, 
and CO Outflow Velocity 

Eigure [14] shows the minor axis position-velocity dia¬ 
gram, with north corresponding to positive offsets. We 
use this to study the minor axis gradient, the outflow 
velocity, and line splitting. Particularly at large vertical 
displacements, the CO is patchy with incomplete cov¬ 
erage, so to build Eigure [14] we integrate over a width 
±1.5 kpc around the minor axis of the galaxy. That is, 
this is the integrated position velocity diagram for the 
outflow region in Eigure [6] 

Purple lines in Eigure [14] show the systemic velocity 
of the galaxy (solid), the mean velocity for all emission 
with | 2 ;| > 1 kpc above or below the disk (dashed), and 
the 95% envelope of emission > 1 kpc above the disk 
(dotted lines). The outflo w is expected to be launched 
from a small central region (|ShoDbell fc Bland-Hawthorn] 
[THH), so that, naively, the difference between the aver¬ 
age velocity above and below the disk and the systemic 
velocity of the galaxy should yield the projected speed of 
the outflow, (^outflow) sin i = (^U| 2 ;|>i kpc) ^sys)* Table 
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Fig. 14.— Integrated CO (2-1) position-velocity cut along M82’s outflow. We plot intensity integrated over ±1.5 kpc along the major 
axis as a function of displacement along the minor axis. Intensity in the image shows intensity of emission with contours stepping by 
factors of two. Green lines show the median velocity for each vertical displacement, revealing increasing velocity with displacement from 
the plane. This signifles either acceleration or a wide range of ejecta velocities. The solid purple line shows the mean velocity while the 
dashed lines show the mean velocity above 1 kpc vertical displacement. The dashed purple lines show the 95% envelope of emission over 
the same range, an indicator of the maximum ejecta velocity. Note the very clear line splitting to the north of the galaxy (discussed in 
terms of a conical geometry, below), but this is less clear to the south, where the lines are simply very broad. 


m summarizes these velocities. The projected mean ve¬ 
locity offset between the high latitude emission and the 
systemic velocity is 78 km s“^. The inclination, i, 
of the outffow is uncertain and because i is small, this 
value exerts a large inffuence on the derived velocity. For 
i = 10°, ^outflow ~ 450 km s“^, while for a moderately 
higher i = 15°, this drops to ^outflow ~ 300 km s“^. Thus 
for our working inclination of i = 10°, the CO outffow 
speed is 450 km s“^. For comparison, M82’s ma ximum 
rotat i on velocity is jus t under 200 km s“^ (see iSofnd 
119981: iGreco et al.l |2()12[) and the Ha outffow velocity is 
525 to 625 km s“^ (|ShoDbell fc Bland-Hawthorn! 1 1998F 
The velocity gradient along the minor axis appears as 
a green line in Figure [HI Southern positions show more 
blueshifted emission while northern positions show more 
redshifted emission, on average. This has the same sense 


as the Ha velocit y gradient (iMcKeith et al.l [l995[) and 
the Hi gradient (|Cottrelll I1977D . which is increasingly 
blueshifted emission below the disk and increasingly red- 
shifted emission above the disk. The magnitude of the 
gradient is ~ 32 km s“^ kpc“^ in the plane of the sky. 

We interpret Figure [14] as a blend of two phenomena, 
rotation along the major axis and outffow along the mi¬ 
nor axis. In a moment, we will see that rotation slows 
with increasing height. Meanwhile, a conical outffow 
morphology produces a double peaked line centered on 
the mean velocity of the outffow, with the two compo- 
ne nts reffecting the different tilts of t he edge (for more, 
see lShopbell fc Bland-Hawthornlll9^ . A blend of these 
two components will produce Figure [Hjand will naturally 
create the gradient. Rotation dominates within ±0.5 kpc 
of the disk, producing the broad, almost contribute most 
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emission. This is also the regime where Figure [TT] show 
a clear excess of emission in the CO and Ha relative to 
starlight. In the north the line splitting is stunningly 
evident in this range, while in the south the dual-peak 
morphology is present from about 2 to 2.5 kpc but then 
only one component persists beyond 2.5 kpc. To the 
south, the cold material traced by dust and gas is mostly 
to the east of the region of hot dust temperatures (see 
Figure ??), so this may indicate only partial covering of 
the outflow cone by high column cold material. 

Figure [14] also plots the envelope of emission, calcu¬ 
lated as the 95^^ percentile of emission sorted by veloc¬ 
ity. This envelope describes emission at a wide range 
of vertical offsets well. The fact that the envelope does 
not curve substantially in or out argues that acceleration 
does not dominate the CO kinematics. That is, even 
though we observe a velocity gradient, we interpret this 
to reflect a changing blend of components rather than an 
acceleration of material as it moves away from the disk. 

9.2. Line Splitting and Conical Geometry 

The lines in Figure [14] are very broad and especially 
in the northern part of the galaxy the position-velocity 
diagram shows multiple components at a given verti¬ 
cal displacement. This can be seen more clearly from 
direct plots of the peak and integrated spectra along 
the minor axis in Figure [15] “Line splitting,” that is 
observing two components separated in velocity along 
each line of sight, has often been sought and invoked as 
a signature of a emission coming from the edges of a 
cone or cyclinder, pyticularly considering Ha (e.g., se e 

iMcKeith et al.l[T9^ IShopbefl fc Bland-Hawthorn]11998[ ). 

In this sketch, emission arises from two sides of a tilted 
cone of emission being expelled from the galaxy. The 
different inclinations of the two sides of the cone yield 
distinct velocities, one at u = Uout sin (i + 0.5 0) and one 
at u = Uout sin (i — 0.5 0), where 0 is the cone opening 
angle. The difference, Au, depends on the inclination, 
the opening angle of the cone and the outflow velocity 
via: 

0 

Av = 2 Uout cosi sin - . (2) 

For the case of M82, cosi ^ 1 and Av « 

2 Uout sin (0.5 0). Even for an outflow that Alls a cone or 
covers the edges patchily, an approximate form of Equa¬ 
tion [2] should still hold. That is, the line width should be 
is proportional to the opening angle and velocity of the 
outflow for a wide range of geometries. 

Eigure [16] accumulates all of the CO and Hi spectra be¬ 
tween 1.5 kpc and 2.5 kpc from the plane in the outflow 
region. That is, it shows the integrated spectra of the 
area that we consider most likely to reflect only cold ma¬ 
terial associated with the outflow. The two-component 
profiles are clearly visible to both the north and the 
south, though the details differ. In the south the double- 
peaked profile is most clearly visible in Hi, while to the 
north the CO shows clearer line splitting. We have al¬ 
ready seen a changing phase as a function of height, here 
we see evidence from the spectral profile that the mate¬ 
rial around the outflow is a mixture of H 2 and Hi with 
the composition varying around the galaxy. 

Eigure [15] shows line widths measured at quarter of 


maximum for peak and integrated spectra along the out¬ 
flow. We report typical high latitude {\z\ > 1.5 kpc) line 
widths for the CO in Table [d] which are ^ 150 km s“^. 
Eor an outflow velocity, of ^ 400 km s“^ and this Au, 
Equation [2] suggests an opening angle ^ 20°. The peaks 
of the two components in the integrated profiles are some¬ 
what closer together, ^ 100 km s“^, which would ar¬ 
gue for 0 13°. These are somewhat smaller open- 

ing angles than those suggested by Ha emission (e.g ., 

iMcKeith et al.l[T9^ : iShopbell fc Bland-Hawthorn]11998 ), 

which range from ^ 30-60°. A lower outflow velocity 
(e.g., due to a higher correct inclination) could remedy 
the discrepancy or our focus on high latitude material 
may And a more columnated outflow. In either case, the 
double peak profiles reinforce that the cold gas at high 
vertical displacement is associated with the wind seen at 
other wavelengths, likely confining the hot wind. 

10. SLOWING ROTATION WITH HEIGHT 

The CO (2-1) line is wider and smoother closer to the 
disk and the velocity held near the disk shows rotation 
(Eigure [2]). This is the other component visible in Eig¬ 
ure [HI a “puffed up” molecular disk around the star- 
burst. If material is ejected from the central, rapidly 
rotating starburst to larger radii, then conservation of 
angular momentum will cause it to slow down. This slow- 
in g rotation with vert ical displacement has been n oted 
bv lWalter et al. (1200211 and lSeaquist fc Clark] (j2 00 1[ ) (see 
also lSofue et al.l[l992[ ). The left panel in Eigure [17] shows 
that this effect is present in our data along the minor axis 
out to ~ 1 kpc, beyond which rotation is no longer visi¬ 
ble. 

Our data have coarse resolut ion compare d to changes 
in the rotation curve of M82 (|Sofuel 119981 ) and so are 
not ideal to measure structure near the disk. Still, we 
can quantify the decline in rotation to some degree. In 
our data, the slope of the best-fit velocity gradient along 
the inner 2 kpc parallel to the major axis drops by as 
much as 100 km s“^ kpc“^ from the midplane to ±1 kpc. 
This is larger than gradients seen in edge-on galaxies 
with exte nded Hi or ionized ga s distributions (such as 
NGC 891 lOosterloo et al.l[200it ). where the magnitude 
of slowing r otation is typically ^ 15 to 25 km s~^ per 
scale height (|Rand et al.ll2008l ). Taking the scale height 
as ^ 0.5 kpc, we And a drop of ^ 30 to 50 km s“^ per 
scale height in M82. 

Eor rotation to slow so quickly with increasing height, 
material must be moving to larger radius very quickly. 
This measurement argues that material is being violently 
flung out of M82 at a wide opening an gle, not only along 
the v isible hot outflow . Eo llowing ISeaquist fc Clarla 
([200l[ ) and I Walter eHI] ([200l ). we also consider a more 
direct version of angular momentum conservation to 
gauge the flow of material in the immediate vicinity of the 
starburst. They suggest that the quantity L = r±v\\ of¬ 
fers a reasonable observational analog of the angular mo¬ 
mentum. Here r± the projected radius along the major 
axis (the axis of rotation) and v\\ the intensity-weighted 
mean velocity along the line of sight. As discussed by 
ISeaquist fc Clark] (|2QQ1[ ) this assumes the observed v to 
map to the tangential velocity and so is is an approxi¬ 
mation. Still the calculation provides useful insight and 
we plot the resulting angular momentum isocontours in 
the right panel of Eigure [T7l 
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Peak Along Minor Axis Integrated Over Major Axis 




LSR Velocity [km s LSR Velocity [km s 

Fig. 15.— CO (black) spectra as a function of vertical displacement from the major axis, {left) The spectrum at the location of peak 
brightness (i.e., an individual line of sight) at each vertical displacement and right integrated over ±1.5 kpc along the major axis. Note 
the shift of the peak, the line splitting at high vertical displacements in the integrated spectrum (mostly to the north), and the significant 
detections of very broad lines even to very large vertical displacements. Light red horizontal lines show quarter-maximum line widths, with 
the average high latitude velocity widths reported in Table [4| Under simple assumptions, the width of the lines combined with an estimate 
of the outfiow speed suggests an opening angle of 20°. In gray, we plot Hi spectra, which have lower quality, for displacements that are 
not contaminated by absorption. 
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Fig. 16.— Integrated and normalized spectra from 1.5-2.5 kpc above and below the disk in CO (thick line) and Hi (filled histogram). 
Both the northern and southern outfiow region show double peaked line profiles, indicating a conical geometry with cold gas bounding the 
hot outfiow. The mixture of molecular and atomic gas in the profile varies with region. In the south the two-peaked shape is more clearly 
visible in Hi but is not clearly visible in the CO. In the north the double peak appears clearly in CO and Hi. Purple lines indicate the 
systemic velocity (solid) and the mean velocity of molecular gas more than 1.5 kpc above and below the disk (dashed). 
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Fig. 17.— {left) Rotation of the disk, traced by mean velocity along the line of sight as a function of offset along the major axis, for 
emission in the midplane and in cuts parallel to the major axis but above and below the midplane. Rotation remains visible but decreases 
in amplitude with increasing height until it is largely imperceptible by 2 ; ~ 1.3 kpc (80"). Note that for edge-on systems like M82, the true 
rotation curve requires careful modeling of the kinematics. Our purpose with this figure is to show that rotation diminishes with vertical 
displacement, (right) Contours of equal angular momentum, L = for regions of bright CO (/ > 10 K km s“^) plotted over the CO 

intensity image. The gray contours show lines along which angular momentum should be approxima tely conserved and so plausible lines 
along which material kicked out of the starburst might fiow (ISeaanist fc Clarl3l200II : [Walter et al.ll2003j . Black arrows illustrate the implied 
wide angle fiow out around the starburst. Expulsion of cold material in many directions helps explain the amorphous appearance of the 
CO and dust maps. 
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Fig. 18.— Mass flow through a plane parallel to the major axis, 
in Mq yr“^, calculated from multiplying the local mass surface 
density of H 2 (from CO), Hi, or Sgas = Sdust/<^DGR W the mean 
outflow velocity. The error bars indicate the effect of changing the 
inclination by ±5°. For comparison, we pl ot the outflow rate cal¬ 
culate d based on observations of hot gas bv IStrickland &: HeckmanI 
(l200^ and a rough estimate of the star formation rate of the sys¬ 
tem. The mass flow rate corresponds to an energy or momentum 
outflow rate via Equations [3] and |4] and the outflow may account 
for ~ 20% of the energy and momentum produced in stars. Tak¬ 
ing dust as a tracer of all cool gas, the mass outflow rate declines 
with an e-folding length of 1.4 kpc; for an outflow velocity of 
450 km s“^ this corresponds to a timescale of 3 Myr. 

The shape of the iso-L curves in Figure [T7| argues that 
the starburst expels cold material over a wide range of 
angles, not only coincident with the hot wind along the 
minor axis. This provides a natural explanation for the 
round, extended appearance of the CO and dust maps 
compared to the columnated hot outflow. The hot out¬ 
flow cavity is very good at “lighting up” associated mate¬ 
rial but material is expelled over a broad range of angles. 

11. FLOW OF MATERIAL, MOMENTUM, AND ENERGY 
ALONG THE MINOR AXIS 

To understand the amount and fate of the molecular 
gas associated with the M82 outflow, we combine what 
we learn from the kinematics and the mass profiles. If 
we consider an outflow moving straight up along the mi¬ 
nor axis and adopt an outflow velocity, then the mass 
profiles in Figure [10] imply a mass flux through each suc¬ 
cessive surface. Figure [THl shows the mass flux through 
progressive surfaces along the minor axis for a fixed out¬ 
flow speed of 450 km s“^. We also experimented with 
using the local velocity field to calculate a mass flux. 
This yields similar results; the approach to calculate the 
velocity appears to matter less than the large system¬ 
atic uncertainty introduced by our imperfect knowledge 
of the inclination. At our adopted speed, gas will travel 
1 kpc in 2.3 Myr. Error bars indicate the effect of chang¬ 
ing the inclination of the outflow by ±5°. Along with the 
uncertainty in the mass estimates already discussed, this 
renders the outflow rate uncertain by a factor of ~ 2- 


3. However, the internal consistency is likely substan¬ 
tially better than this. For comparison, the Figure also 
shows the approxima te star formation rate of the galaxy 
(|Westmoqnettell2QQ'^ and the mass outflow rate in ferred 
for the hot wind by iStrickland fc HeckmanI (j2QQ9[ ) . 

FigureUHlcan be recast in terms as a rate of momentum 
or energy flow away from the disk. For a fixed outflow 
velocity, the corresponding energy rate is 


dE ^,^40 _i / M 

—— = 2x10^^ erg yr ^ - t 

dt I 10 Mq yr ^ 


and the momentum flux is 


450 km s ^ 
(3) 


^ = 4500 M© km s Ur ^ 
dt \ 10 Mq yr 


-1 


450 km s 


-1 


(4) 

Figure [Island its momentum and energy analogs show 
first that cold material seems to be leaving the outflow as 
it moves away from the disk. That is, the profiles decline. 
Second, the outflowing cold gas is changing phase to be¬ 
come more atomic and less molecular; the CO (blue) 
profile declines fastest. Third, the amount of cold gas 
near the disk is so large that the inferred rates dwarf the 
hot gas outflow rates near the disk. However, the outflow 
rate of molecular gas, at least, is less than or comparable 
to the hot gas outflow rate by « 4 kpc. Thus despite the 
visible molecular outflow and large amount of CO emis¬ 
sion, the molecular wind is not a major factor on large 
scales. Fourth, mass flux rates around the disk are of the 
same order or even larger than the star formation rate 
in the galaxy. Thus the redistribution of material due to 
stellar feedback will be just as important to the evolu¬ 
tion of the starburst as the depletion of gas into stars. 
Finally, the figure shows that distinguishing outflowing 
from tidal Hi and dust at large distances is critical to 
understandi ng the mass budge t in th e outflow. 

Following lOstriker fc Shett^ (|2Q11[ ) , a reasonable mo¬ 
mentum output due to supernova is « 3000 Mq km s~^ 
per 1 Mq formed and a reasonable kinetic energy in¬ 
jection is ~ 10^^ erg per 1 Mq formed. For a star 
formation rate of 10 Mq yr“^ there should then be 
^ 30,000 Mq km s~^ yr“^ and 10^^ erg yr“^ released 
into the ISM. Above about 1.5 kpc (where we are fairly 
confident that the outflow dominates), the momentum 
associated with the mass outflow rate of 10-20 Mq yr“^ 
represents ~ 10-30 % of the momentum available from 
star formation each year and the energy available repre¬ 
sents a similar fraction of the available energy. Both the 
energy and momentum injection and the outflow num¬ 
bers are uncertain, but this first order comparison shows 
plausible agreement. 

Again, we note that this translation of distance from 
the disk into evolution of the outflow with time assumes 
a roughly constant mass outflow rate over the last ^ 
10 Myr. A useful future avenue of investigation, but one 
beyond the scope of this paper, would be to attempt 
to combine these measurements with information about 
the recent star formation history in M82’s bu rst (e.g., see 
Ide Griisl I2QQ1I: iForster Schreiber et al.l |2QQ3[) along with 
plausible lag times for launching of a wind. 
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12. SYNTHESIS 

We present the most sensitive wide-field map of the CO 
emission around M82 to date, which traces the distribu¬ 
tion and kinematics of molecular gas in the sprawling 
gaseous halo that surrounds the famous starburst. This 
includes extended CO emission coincident with other 
tracers of M82’s well-known galactic wind, with clear de¬ 
tections at distances > 3 kpc above the starburst disk 
and well outside the stellar disk. By co mbining these 
new CO data with maps of du s t emission (iRoussel et al.l 
I 2 MIIE ngelbracht et al.l I2QQ6I: iKaneda et al ]M), Hi 
(lYun et ^ I 1993 I 119941 ) , tracers of the ho t out¬ 


flow (iLee et al.l I2QQ9I: iStrickland fc HeckmanI [200^ 
iHoopes et al.ll2QQ5f) . starlight (jJarrett et a .11200311. and 


other CO transitions (C O 1-0 and CO 3-2 1 Walter et al.l 
l2QQ2l : IWilson et al.ll2Q12l ) we attempt to build a complete 
picture of the cool gas and dust halo around M82. 

Figure [ 19 ] illustrates our conclusions schematically. We 
And CO associated with the hot outflow. Both CO and 
Hi play an important role confining the outflow. The in¬ 
terface between the hot wind cavity, with its associated 
intense radiation held, and the surrounding cold mate¬ 
rial gives rise to many of the most striking aspects of 
the outflow: the Ha and PAH emission, the hot dust, 
and scattered FUV light. The association of the cool 
gas with the outflow is clearest outside ^ ±1.5 kpc from 
the disk. In these extended regions double peaked line 
profiles indicative of conical geometry and a minor axis 
velocity gradient both link the gas to the outflow kine¬ 
matically. The coincidence of CO, hot dust, and Ha 
along the minor axis also support this picture. 

However, although cold material is associated with the 
hot outflow, it does not appear to travel far into the halo. 
Mass profiles of dust, CO, and Hi all decline along the 
minor axis without any corresponding evidence for accel¬ 
eration. Thus basic mass continuity arguments suggest 
that cold material is leaving the outflow, stalling out or 
falling back to the disk. Dust should remain mixed with 
any low temperature ionized gas, so material vanishing 
into an ionized phase does not appear to offer a plausible 
alternative. There is, however, a clear phase change as¬ 
sociated with the molecular gas, which tends to become 
more atomic as it flows out of the galaxy. 

A wide held view of dust, gas, and stars around M82 
supports the idea of a cold “fountain” rather than a cold 
outflow. Profiles of dust emission over a wide field do 
not diverge. That is, their profiles appear significantly 
steeper than a projected p oc profile, which would be 
required for mass continuity and a fixed outflow veloc¬ 
ity. IR colors suggest that a changing dust temperature 
does not offer a plausible alternative scenario, the dust 
seems to stay warm enough for HerscheFs SPIRE instru¬ 
ment to detect it to large distances. Most dust, and thus 
presumably most of the cold material ejected from the 
galaxy, is concentrated within a few kpc of the galaxy. 
Its morphology at low intensity is approximately circu¬ 
lar in projection, suggesting material is not just ejected 
along the hot outflow cavity. 

Expulsion of material in all directions is supported by 
the CO kinematics closer to the disk. Slowing rotation 
with height and a direct examination of lines of fixed 
angular momentum show “flow lines” with much wider 
opening angles than the hot outflow. This is consistent 


with cold material being ejected in many directions but 
preferentially lit up along the hot wind cavity. 

Thus we argue that the data support a highly dynamic 
cold gas fountain around M82. Material is hurled out as 
molecular gas but becomes atomic over a timescale of a 
few Myr. At the same time, material stalls out or falls 
back on a slightly longer timescale. Doing so it must 
change phase and position around the galaxy, because 
fallback is not readily evident from the CO minor axis 
position-velocity diagram. 

The mass flux along the minor axis changes accord¬ 
ingly. It is much higher than the outflow rate associ¬ 
ated with the hot wind or the star formation rate near 
the galaxy, but it drops with increasing distance. By 
« 4 kpc from the disk, the mass flux of molecular ma¬ 
terial is not an important contribution compared to the 
hot wind. The situation with the total or atomic gas is 
less clear. The mass flux also drops, but in future work it 
will be important to distinguish tidal features from out¬ 
flow features in the atomic gas to cleanly measure the 
outflow component. Still, the implied trend, especially 
considering the dust to be ejected material, is that the 
cold gas outflow does not reach far out into the halo, un¬ 
like the hot material. The motion of this cold material 
may account for 10-30% of the energy and momentum 
produced by star formation each year. 

These results sound something of a cautionary note for 
outflows detected in unresolved or marginally resolved 
spectra. We show in the appendix that if viewed at a dif¬ 
ferent inclination, M82 could produce a dramatic exam¬ 
ple of the wide CO line wings commonly used to identify 
outflows. However the evidence we present argues that 
the molecular material, in particular, is not a true out¬ 
flow. Large scale feedback effects on the molecular gas 
are clearly critical to the evolution of the system. They 
will remove gas from the starburst region and affect the 
structure of the galaxy but we do not see good evidence 
that molecular gas makes it far out into the halo of M82. 

Einally, we note several areas for future investigation 
that will help shed further light on this key example of 
a molecular wind. We have assumed a constant recent 
mass outffow rate, but in fact it will be productive to 
jointly model the ffow of material out of the galaxy and 
the recent star formation history of the galaxy. Also, 
although we present good arguments that radiative and 
not collisional heating dominates the dust energy bud¬ 
get, the role of shocks and hot gas heating dust at the 
interface of the hot superwind and cold material has yet 
to be explored in detail. Last, the edge-on orientation 
that makes M82 so tractable to this sort of study also 
leaves the inclination of the galaxy uncertain; any more 
stringent constraints on this quantity will have impor¬ 
tant implications for how we think about the energetics, 
momentum, and possible escape of the material leaving 
the burst. 


We thank the anonymous referee for a careful and con¬ 
structive report that improved this paper. This work 
benefited from helpful discussions and feedback from 
many colleagues including Todd Thompson, Eve Os- 
triker, Jeffrey Mangum, John Hibbard, and the Uni¬ 
versity of Virginia/NRAO star formation group (Kelsey 
Johnson, Amanda Kepley, and Scott Schnee). This pa- 
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Fig. 19.— A schematic of the view of cool gas and dust around M82, not to scale. Blue notes cite evidence for the red statements. 
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APPENDIX 

PROSPECTS FOR DETECTING OUTFLOWS SIMILAR TO M82 IN UNRESOLVED SPECTRA 

In both the radio and the optical, unresolved spectra remain a main tool for identifying candidate outflows. Though 
it may be the best-known galactic wind, the M82 outflow would be difficult to identify from integrated spectroscopy 
because of its high inclination. As Figure [2Ql shows, projection effects render the wind almost invisible in the global 
spectrum of the galaxy. By measuring individual spectra for our sub-regions (Figure [6]), we simulate the effects of 
seeing M82 from other angles. To do this, we assume the outflow to be orthogonal to the disk and that the true 
inclination of the galaxy (outflow) is i = 80° (10°). For a series of other possible inclination angles, igim, we adjust 
the width of the disk spectrum by a factor of sini/ sinigim while conserving its integrated flux. We carry out a similar 
operation on the outflow and add the two to construct a new, artificially inclined spectrum. Mechanically, we achieve 
the adjustment for inclination by simply interpolating the old spectrum, with the abscissa rescaled, on to our fiducial 
velocity grid and rescaling the interpolated spectrum to match the flux of the original version. We expect this approach 
to fail at inclinations below igim ~ 20° because the intrinsic velocity dispersion within the starburst, which we do not 
model, will become important. Therefore we consider igim over the range 20 to 80°. 

The top left panel in Figure [20] shows the integrated spectrum of M82 as we observe it. The top right panel shows 
the integrated spectrum simulated for different inclinations. The spectral s ignature of the outfl ow in a face-on version 
of M82 is a broad, low component (similar, e.g., to that in Markarian 231. iFernglio et al']l2010[ ). We quantify the flux 
and amplitude of these faint line wings in the bottom panels of Figure |20l For these calcul ations, we define t he “line 
wings” to be all channels that have flux density less than 5% of the peak value (similar to lLerov et al.ll^Qlbf ). For a 
Gaussian line profile, just over 1% of the flux emerges from the line wings defined in this manner. For the flux in 
these wings, we sum all channels that meet this criteria. To measure the average amplitude, we calculate j ^Fv 
where F^ < 0.05 that is, we derive the flux-weighted flux density in part of the spectrum with amplitude below 

5% of the peak value. 

Figure 120] shows results for two cases: (1) the entire outflow region and (2) a spectrum constructed by doubling the 
emission outside ±1.5 kpc. At these high latitudes, we expect to avoid foreshortened disk material and gas ejected at 
a range of angles, observing only material flowing away from the disk along the minor axis. We double the spectrum 
because from mass continuity arguments (see main text), we expect that at least an approximately equal amount of 
material must flow out along the minor axis inside ±1.5 kpc. Therefore, doubling the spectrum represents a realistic. 
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Fig. 20.— Prospects for detecting M82-like outflows from unresolved CO spectroscopy, (top left) The integrated spectrum of M82 (black) 
from our masked data cube with the spectra of the individual regions — the starburst disk (red), the whole outflow (purple), the outflow 
outside z ± 1.5 kpc multiplied by two (blue; see text), and the extended disk (gray). Conserving flux and applying varying corrections for 
inclination, we estimate how the integrated spectrum would appear if we viewed M82 from other angles {top right). The bottom left panel 
shows the fraction of flux in the “wings” of the integrated spectrum, with wings deflned as channels with < 5% of the peak flux density. 
We plot the values for a Gaussian line profile and the real integrated spectrum as horizontal lines for comparison. The bottom right panel 
shows the amplitude of the line wings in the simulated spectrum; here amplitude is deflned as the flux weighted flux density in parts of the 
spectrum with flux density < 5% of the peak value. Results for only the high latitude material suggest a conservative lower limit of 5% 
flux in line wings with magnitude ~ 1% of the peak flux density for a moderately inclined version of M82. 


but still conservative (because it accounts for no phase change or fail-back) lower limit for this exercise.This doubled 
|z| > ±1.5 kpc curve represents our conservative lower limit and the main result of this exercise: in a face-on version 
of M82, faint line wings with amplitude < 5 % of the peak are likely to contain at least ~ 5% percent of the total 

line flux. The results for the whole outflow region suggest that the true signature of a face-on M82 could be even 
stronger than this and could be as high as > 10-20% of the flux in line wings with flux density ^ 3% the maximum. 
These more optimistic results carry the caveat that they may include blended disk emission and emission from ejecta 
at a wide range of inclinations. For these ejecta a a wide range of angles, the effect will be to move some of the flux 
to a narrower component with higher amplitudes, one that may be harder to pick out from the disk in an integrated 
spectrum. 


THE CO-TO-H 2 CONVERSION EACTOR IN THE OUTELOW 

The procedure to convert our observed CO luminosity into a molecular gas mass is complicated by the fact that 
physical conditions in the o utflow may differ substantia lly from those in the disk galaxy giant molecular clouds usually 
used to calibrate aco (e.g.. I Donovan Mever et al.ll2Q13[ ). In this section, we step through what we can infer about the 
conversion factor in th e M82 outflow based o n comparing different CO transitions and tracers of different ISM phases. 
We refer the reader to iBolatto et al.l (j2Q13b[ ) for a thorough discussion of all of the issues involved in estimating the 
CO-to-H 2 conversion factor. 


Optically Thin Limit 

Following the treatment of CO ( 1 - 0 ) in the NGC 253 wind bv IBolatto et al.l (j2Q13a[ ). we consider the case of an 
optically thin line in local thermodynamic equilibriu m (LTE) as a lower limit. Using Einstein A coefficients and other 
data from the Leiden Atomic Molecular Database (jSchoier et al.l [ 200 ^ . we calculate the optically thin aco for the 
three lowest J CO transitions: 
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Her e we have assumed local thermodynamic equilibrium (LTE) and neglected CMB effects; see iBolatto et al.l (|2013aD 
and iLerov et al.l (|2014l) for an explicit list of the equations used in the LTE calculation. The implied line ratios in 
the optically thin case, R 21 = hi/ho ~ 2.8 and i ^32 = h 2 /h^ ^ 1- 3 (both in brightness temperature units) for 
Tex = 30 K match expectations from RADEX (|van der Tak et ahl 120071) in the limit of large line width, low column, 
and high volume density with Tkin = 30 K. Eor the LTE calculation, varying Tex from 20-40 K makes only a ^ 20% 
difference on The impact is larger for a non-LTE treatment, as we see below. 


Comparing Low-J CO Transitions 

The line ratios implied by Equation IB II are high, with both R 21 and T 32 greater than unity; indeed R 21 ^ 3. M82 
has been studied extensively in molecular gas, so t hat we have some knowledge of the true line ratios. Across the 
whole system, comparing our luminosity to that of iSalak et al.l (|2013[ ) yields R 21 ~ 0.8. Though the apertures are 
imperfectly matched, dividing our outffow luminosity by their high-latitude luminosity yield s T 21 0>3 but cert ainly 

not above unity. M ore concretely, we have in-hand CO ( 1 - 0 ) and CO (3-2) maps from IWalter et al.l (| 2002 f ) and 
IWilson et al.l (| 2012 [) . We use these for a direct calculation of line ratios, although note that this comparison is limited 
by the coverage and sensitivity of the other data (see Eigure [ 2 T]) . 

To compare our CO ( 2 - 1 ) map with the other transitions, place all three CO data cubes on to a matched grid with 
coarse ~ 40" x 40" x 20 km s“^ pixels. The centers of the independent pixels appear as circles in Eigure [21] along with 
the coverage of the CO ( 1 - 0 ) and CO (3-2) maps. In Eigure [211 we plot the spectra for transition at each pixel with 
intensities (in Kelvin) normalized to the peak of the CO (1-0) spectrum. We also plot the two line ratios R 21 and 
T 32 , color coded by the 250/im intensity in the pixel, which serves as a rough indicator of distance from the starburst. 

These spectra show that the overall line profile shape remains the same among the transitions, although the signal 
to noise of the other data are lower than our CO ( 2 - 1 ) map (this is particularly an issue before the heavy rebinning 
used to make the figure). This consistency among all three maps offer some check on the accuracy of our calibration 
and error beam correction scheme. In M82, the intensities of the CO ( 2 - 1 ) and (3-2) lines tend to be lower than 
that of the CO ( 1 - 0 ) line, but not by a large factor. The exceptions to this general agreement are the corners of the 
grid, w here the fini t e exte nt of the other maps make the comparison unreliable. Euture comparison to the wide field 
map of ISalak et al.l (|2Q13D will allow a check on whether the variations reflect limitations of the data or v ariations in 
physic al conditions in regions of faint, extended emission. Line ratios near unity but not above agree with lWeifi et al.l 
(l2QQ5f). who s tudie d multiple CO transitions for individual pointings across the streamer/outflow area as defined in 
I Walter et al.l (j 2 QQ 2 [ ) and found average R 21 = 1.0 ± 0.1 and T 31 = 0.8 ± 0 . 2 . Eor comparison, we find median R 21 = 0.7 
with « 0.15 scatter and median T 32 ~ 0.75 with 0.3 scatter. Eigure!^ shows that the ratios are reasonably indicative 
of areas across the overlap between the maps, though the central burst may be slightly more excited. 

These ratios are not what one would expect for optically thin gas in LTE. To investigate what physical condi- 
tions could produce the observed ratios, we compare these ratios to a large grid of Radex calculations. Radex 
(|van der Tak et aP I2007D is a non-LTE code that models line emission using an escape probability formalism for 
a single set of physical conditions (density, temperature, and column density per line width). We calculate expected 
^^CO emission for a large grid of these physical parameters and compare our observed line ratios to these calculations 
in Eigure [ 22 l Our interest here is broad results, not a point-by-point comparison, so we identify points in the grid 
that match our median T 32 and R 21 with a tolerance of ±0.2. This does a good job of capturing most of the observed 
ratios (Eigure [ 2 T]). Eigure [ 22 ] shows our two main results for the calculation. Eirst, almost all models that match 
the observed ratios are optically thick in all three low-J CO lines, further arguing against the optically thin case. 
Second, the observed line ratios can be produced by either cold (Tkin ^ 20 K) and dense (log^^g ^H 2 [cm~^] > 3.5 gas 
or hot (Tkin ^ 30 K) but rarer (log^^g ^H 2 [cm~^] ^ 2 gas. These results agree with IWeifi et al.l (j2005[ ). who found 
significant opacity in the lines considering individual pointings in the same area we study here. By including 

t hey were ab l e to fu rther constrain log^^g ^H 2 [cm~^] ^ 3 cm“^, towards the lower end of viable densities. This 
also led IWeifi et al 1 (l 2 nn^ to argue for a higher tem perature, Tkin 30 K. Note, also that this temperature cannot be 
ar bitrarily high: in ad dition to the constraints from IWeifi et ()2005[ ) , direct observations of the H 2 rotational lines 

by iBeirao et al.l (|2015[ ) constrain the bulk H 2 temperature to be less than a few hundred K and even this “cooler” 
component of warm H 2 represents only a modest fraction of the total mole cular gas niass in the system. 

What does opaque, warm gas mean for the CO-to-H 2 conversion factor? IWeifi et al ] (IWili used analysis of a more 
limited set of CO lines, but still including optically thin tracers, to demonstrate a reasonable correspondence of aco 
with the expected for an opaque, virialized cloud. Eollowing this logic, a Tkin ~ 2-3 times higher than Milky 
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Fig. 21.— Comparison of different CO transit ions in M82. (top left) SPIRE 250 am map with conto urs indicating the coverage of the CO 
(1-0) OVRO+FCRAO map (| Walter et al.ll200'^ ') and the JCMT CO (3-2) map (| Wilson et al.ll201^h Dots indicate the central positions 
for the heavily binned spectra in the bottom panel, (top right) The three line ratios as a function of SPIRE intensity, which roughly 
corresponds to displacement from the central starburst. Here each point is an individual line of sight at the matched 20^^ resolution of the 
HERA map. (bottom) Spectra for all three low-J CO transitions binned to 40'' independent pixels and 20 km s“^ channel width to 
increase the signal to noise. The grid cells are centered at the dots indicated in the top left panel and the average 250/rm intensity in each 
bin is noted in that panel. 


Way clouds, with a similar density (see iHever et al.ll2QQ9D implies a conversion factor ^ 1/3-1/2 the “disk galaxy” 
value. This roughly matches the values found bv lWeifi et al.l (|2QQlf ) in positions that were vertically offset from the 
starburst disk, « 5-10 x 10^^ cm“^ (K km or « 0.25-0.5 times the Milky Way value. It also matches the value 

of ~ 10^^ cm“^ (K km s“^)“^shown bv iKeto et al.l (|2005l ) to produce approximate virial equilibrium in the resolved 
cloud population in the M82 starburst. 
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Fig. 22.— Results for RADEX modeling of line ratios around M82. We run a large grid of RADEX calculations covering a 

range of kinetic temperature (Tkin), collider density (nii2)i and column per unit line width (Nqo /These two panels show our main 
results: (left) the observed line ratios are consistent with models that have either low temperatures (Tkin ;$ 30 K) and high densities 
(log^^Q nii 2 > 3 cm“^) or high temperatures (Tkin ^ 30 K) and low densities (log^^Q nH 2 ~ 2-3 cm“^). Red and blue lines show isocontours 
of the CO (3-2) / (2-1) (blue) and (2-1) / (1-0) (red) ratios for a Ac-o/Au 10^^ cm“^ (km s“^)“^, a reasonable value for a molecular 
cloud. Gray points show all combinations of T and n in our grid that match the typical observed line ratios within 0.25 (we add a small 
amount of noise to the data for plotting), {right) Optical depth for all grid solutions that match the observed line ratios. The calculations 
suggest optica lly thick emission in all three lines, consistent with the comparison to LTE. Work combining ^^CO lines with optically thin 
transitions bv lWeifi et al.l (120051) suggests the high Tkin, low nH2, optically thick solutions as most likely, which in turn suggest conversion 
factors ~ 2-4 times lower than Galactic. 



logto I 350 [MJy sr ‘] 



Fig. 23.— CO (2-1) emission (/co) as a function of far-infrared intensity at 350/rm (/350) measured with SPIRE on Herschelior a reference 
sample compiled from the HERACLES and KINGFISH surveys (contour and histogram, taking all lines of sight with /350 > 10 MJy sr“^) 
and the three regions that we use to analyze M82. Though the effects of Hi and metallicity are clearly visible at low intensities and as a 
tail in the histogram, the reference sample suggests a typical ratio of log^^o-^co/-^350 = —0.75 ± 0.2 for the molecule-dominated parts of 
disk galaxies. M82 shows a systematically higher amount of CO emission per unit IR emission, log^^Q/co/-^350 = —0.35. The translation 
of /350 to dust mass is not trivial, but the simplest explanation for the difference is a conversion factor, aco that is 2.5 times smaller in 
M82 than a typical disk galaxy. 


Comparing Dust and Gas 

Dust emission offers an powerful, orthogonal constraint on the conversion factor. Following [Israeli (jl997l ). the dust 
mass derived from the infrared spectral energy distribution has been used as a tracer of the total neutral gas column. 
This technique has been used to estimate H 2 mass without the use of CO in a variety of sy stems. Indeed, versions 
of this approach have been appli ed to the inner region o f M82 as early as IWild et ahl (jl992[ ). Details of using dust 
to solve for aco are discussed in iSandstrom et aTI (j2Q13f ). iLerov et al.l (j2Qllf ). and references therein. Here we first 
compare the correlation between CO and long wavelength emission measured by Herschel in the outflow region to 
that observed across a wide set of nearby galaxies. We demonstrate that CO emission is bright compared to 250/im 
emission in M82, a strong empirical indication of a different aco- Then we compare the dust surface density derived 
from SED modeling, Edust, to Hi and CO emission to estimate aco- We solve for the dust-to-gas ratio, (^dgr, in the 
extended molecular disk and then, assuming this value to be approximately constant, apply it to constrain gco in the 
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outflow region. 


Suhmillimeter Intensity and CO Emission 

Figure [23] /co as a function of 350/ini intensity measured by Herschel^^. The 350/im intensity, though not entirely 
on the Rayleigh-Jeans tail, will be approximately proport i onal to Td„st x S dnst and for a constant dust-to-gas ratio, 
/350 oc Tdust X ^gas (foi* caveats see iSandstrom et ah! 120131: iLerov et~al]l 2 Qll[) . Similarly, in the limit where molecular 
gas dominates the gas reservoir, Sgas ^ <aco x /co, so that the ratio Ico/hbo is an observational “color” closely 
related to the CO-to-H 2 conversion factor. If all properties other than aco remain fixed, then ac o oc (Ico/hr^n)~^ i n 
regions where H 2 dominates the gas reservoir (see related disc ussion cast in terra s of luminosity in [Groves et al.ll2Q15f ). 

As a benchmark, we t ake th e entire HERACLES sample fsee [Lerov et al.[r2Q13[) and the Herschel image s from KING- 
EISH [Kennicutt et al.[ (|2Q11[ 1. We convolve HERACLES to the resolution of the 350/im data ffollowing [Aniano et al.[ 
Imii) and extract matched intensities for each line of sight where /350 > 10 MJy sr Eigure [23] shows the re¬ 
sulting correlation appears as a data density plot in /co -^350 space and a histogram. The figure also motivates the 
adopted intensity cutoff; below /350 MJy sr“^ the effects of metallicity and Hi on the ratio /co /^350 become sub- 
stantial. That is, the d ust-to-gas ratio begins changing and Hi constitutes a non-negligible part of the gas reservoir 
(|Sandstrom et al.[[2Q13[) . Indeed, the Eigure shows that these effects are already visible, though not yet dominant, 
below /350 ^ 30 MJy sr“^. Erom either the ridge line (median Iqo at a given / 350 ) in the left panel or the histogram 
in the right panel, we calculate the ratio for bright emission in the HERACLES and KINGEISH overlap to be 

logio = —0.75 ± 0.2 where /350 > 30 MJy sr“^ (B5) 

with 0.2 dex the rms scatter among the data, not the uncertainty in the mean and the relationship calculated for 
mostly star-forming disk galaxies at ^ 1-2 kpc resolution. 

Eigure [23] also plots our measurements of M82. We find a good correlation between /350 and /co in M82 but also that 
the normalization of the M82 correlation differs from the reference sample. Outside the highest and lowest intensity 
regions, M82 exhibits systematically higher /co /^350 than normal disk galaxies. This is exactly the sense expected for 
a low aco- Ignoring for the moment the low and high intensity deviations, the median log;L 0 ^co2-i/735o MS2 is 
—0.35, which is 0.4 dex or a factor of 2.5 times lower than the median for the reference sample. 

Taken at face value, /co2-i/735o implies that in M82, is ~ 2.5 times lower than in a typical disk galaxy. 

Effects other than aco can alter this ratio, however. The main effects should temperature variations (/350 oc T^ust, the 
contribution of Hi to the total gas reservoir, and potential differences in the dust-to-gas ratio between M82 and the 
reference sample. The sense of a high Tdust in M82 would be to make /350 higher, which would only lower the inferred 
cteo- Similarly, any additional gas component would render this simple estimate an upper limit. Only the dust-to-gas 
ratio represents a plausible way to recover cteo nearer to the overall reference sample. However, the weaker shie lding 
associated lower dust-to-gas ratio also typically suppresses CO emission (e.g.. [Lerov et al.][2Qlll: [Lee et al.[[2Q15L and 
references therein), which would act in the direction of again lowering /co2-i/735o- 

CO ( 2 - 1 ) emission thus appears about 2.5 times brighter relative to 350/im emission in M82 than in the whole 
HERACLES sample. Though this is an empirical indicator with caveats, we argue that the most straightforward 
interpretation of this high ratio is a lower CO-to-H 2 conversion factor in M82 compared to most of the bright regions 
in HERACLES. At very low and very high /350 this conclusion appears less secure, likely due to the importance of Hi 
at the low intensity end and the complexity of emission from immediate starburst region (which is not the topic of 
this paper) at the high intensity end. 

Results Using Dust Speetral Energy Distribution Eitting 

Eollowing [Sandstrom et ah] (j2Q13[ ) and [Lerov et al.[ (j2Qll[ ). we adopt the more physical — but also more model- 
dependent — approach of comparing the Edust inferred from t heoretical models t o Iqo and Ehi- As described in the 
main text, after cleaning the infrared maps we fit the models of [Draine fc Li] (j2QQ7[ ) to the infrared SED point-by-point 
at the resolution of the SPIRE 350/im map. We compile matched vectors of /co, Ehi, and Edust at this resolution 
and attempt to s olve for aco by minimiz ing the scatter in the dust-to-gas ratio as a function of aco- 

As discussed in [Sandstrom et HI (j2Q13[ ). the method requires both Hi and H 2 to contribute meaningfully to the ISM 
and a single dust-to-gas ratio to describe the system. In M82, this condition is best met in the extended disk and 
tidal streamers (the gray region in Eigure [6]). Eigure [24] shows the dust-to-gas ratio as a function of tested aco in this 
region. Eor several sub set of the data for this re gion, we obtain ^ 2.5 Mq pc“^ (K km s“^)“^, just over half the 
average value found bv [Sandstrom et ^ (j2QI3[ ) using the same CO line and the same dust models for a large sample 
of disk galaxies, = 4.4 Mq pc“^ (K km s“^)“^. 

Thus in the extended disk of M82, scatter minimization in the dust-to-gas ratio implies aco about half the value 
found in normal disk galaxies. This this is approximately the same result returned by comparing infrared intensity at 
350/im to CO emission and so reinforces that simple result. More, the fit yields a dust-to-gas ratio for M82 that we 
expect to also hold in the cold gas entrained in the wind and likely the starburst. The dust-to-gas ratio implied by 
our best-fit aco is Jdgr ~ 0.009, or about IIO times more gas than dust by mass. Eor comparison. [Sandstrom et ah] 


Note that the results are essentially the same if we instead 
use 250/rm emission; we prefer /350 only because it lies closer to 


the Rayleigh-Jeans regime. 
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Fig. 24.— Constraints on aco from comparison of CO emission to dust and Hi. (top left) Scatter in the dust-to- ^as ratio as a function o f 
our adopted aco for th e extended disk, where we have sufficient dynamic range in H 2 /H 1 to solve for aco foUowing iSandstrom et aP (|2013l ) 
and lLerov et al.r(|201in . Our solution in the extended disk is stable across several regimes in CO brightness (different colors). The implied 
dust-to-gas ratio (top right) using our best prescription (black histogram) shows less scatt er and is more consisten t with a single dust-to-gas 
ratio than if we applied the typical DGR from disk galaxies using the same dust models (|Sandstrom et al.ll2013lb The implied dust-to-gas 
ratio in M82 is somewhat lower than in more massive disk galaxies, perhaps consistent with its lower mass. Bootstrapping this dust-to-gas 
ratio to the rest of the galaxy, we constrain in the area around the starburst and the outflow. The resulting scaling between dust and 

total neutral gas (H 1 +H 2 ) appears in the (bottom right) panel, which uses a two-component prescription for aco (Equation [BS]). Outside 
of the brightest part of the starburst, th is prescription yie l ds a r easonably fixed dust-to-gas ratio across the system, achieving the same sort 
of linear scaling of dust and gas seen bv ISandstrom et al.l (|2013l l for a large set of disk galaxies (red points and line). If a single dust-to-gas 
ratio is expected for M82, this approach is a clear improvement over applying only a typical aco (bottom left), which yields larger scatter 
and a non-linear relation between gas and dust. 


(|2Q13[ ) found \ogiQ Sjjgr = —1-86, or l-to-72, for a large set of disk galaxies using the same line and models; that 
study focused on the CO-bright inner parts of galaxies, so this average value may be so mewhat high c ompy ed to, 
say, the Solar Neighborhood or integrated galaxies. Considering integrated galaxy SEDs, iDraine et al.l (|2QQ7[ ) found 
logiQ 6dgr ~ —2, or 1-to-lOO using CO (1-0) and the same dust models. Thus the extended disk of M82 shows a 
slightly lower dust-to-gas ratio than metal rich disk galaxies, but well within the range of normal values seen in previous 
studies using similar data and perhaps consistent with M82’s low stellar mass. We emphasize com parative statement 
because several lines of evidence suggest that the normalization of the dust masses produ ced by the lPraine fc Lil (|2QQ7[ ) 
models may be as much as a factor of ^ 2 times too high (e.g., iFanciullo et al.ll2Q15f ). Thus our best estimate for 
the true value of the dust-to-gas ratio in M82 would be ^ l-t o-22Q, slightly lowe r than Sdgr inferred in the Solar 
Neighborhood from depletions l-to-150, e.g., see synthesis in iDesert et aLlfloM) . 

Applying the best first ^dgr from the extended disk, we can place constraints on the molecular mass, and so 
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TABLE 5 

Results for from Dust Comparison 


Region 

for (^dgr — 0.009 

for (5 dgr = 0-01 

for (5 dgr = 0.009, no Hi 
(upper limit) 

based on Ico/hbo^ 

Extended disk 

Outflow 

2.7 

2.2 

< 6.2 

2.5 

... whole region 

1.3 

1.0 

< 2.6 

1.6 

... faint (/co < 5 K km s“^) 

3.4 

2.7 

< 8.2 

4.4 

... bright (/co > 5 K km s“^) 
Starburst Disk 

0.9 

0.8 

< 1.8 

1.1 

... whole region^ 

3.2 

2.9 

< 3.6 

3.2 

... faint (/co < 50 K km s“^) 

1.4 

1.2 

< 2.3 

1.3 

... bright (/co > 50 K km s“^)^ 

3.4 

3.1 

< 3.7 

3.5 


Note. — Units of Mq pc ^ (K km s ^throughout. 

^ Assuming that = 4.4 Mq pc“^ in the reference sample following ISandstrom et~a D dinii and scaling by the difference 

between /coZ-^sso in the region and the reference sample. 

^ We expect that in the bright part of the starburst disk our modeling fails. But the general result of /co2-i/-f350 similar to a 
disk galaxy holds. 


in the outflow and disk. We calculate aco itt three ways: 


1. We take ^dgr = 0.009 and estimate ctco via the equation: 

X] ^DGR X ^dust 


CtGO = 




(B6) 


with the sum over regions of interest: the outflow region, the starburst disk, and the extended disk. We also 
subdivide those regions to explore ctco variations within them. 


2. We also place a conservative limit on aco by neglecting any contribution from Hi: 


ggo < 


Z ^DGR X Ddust 

Z^co 


(B7) 


which avoids uncertainty about the relative geometry of Hi and CO and does not rely on the uncertain Shi 
closer to the starburst disk, where the line goes into absorption. 


3. We apply Equation IB6I but with ^dgr 
dust-to-gas ratio. 


0.01 (|Draine et al.l [200^ as a check on the impact of the adopted 


Table [5] summarizes these calculations. Neglecting the very bright center of the galaxy we distill the findings as 
follows: dust emission suggests that the conversion factor throughout M82 is likely about half the disk galaxy value in 
the extended disk and faint parts of the wind. It appears to be lower, by about an additional factor of 2, in the bright 
regions near the starburst. In the very brightest parts of the starburst itself the ratio /go 2 -i /^350 actually resembles 
that found in the disks of normal galaxies; however, this region is not the subject of the paper. In agreement with the 
multiline comparison, even the lowest plausible gco values from the dust treatment are still many times the optically 
thin value. 


Synthesis 

All of the arguments considered lead us to expect a moderately low conversion factor throughout M82. As a 
co nservative single number , we take the ~ 2 Mq pc“^ (K km or about 1/2 the disk galaxy value argued 

bv ISandstrom et al.l (j2QI3f ) and about 1/3 the value adopted by other papers studying HERA maps of galaxies (e.g., 
iLerov et~ani2QI3[) . This single value agrees within a factor of ^ 2 with almost all resul ts: line ratios, co mparison to 
sub-mm intensity, comparison to dust model, and cloud virial masses in the starburst (|Keto et al.ll2QQ5[ ). It implies 
optically thick emission, which is consistent with the observed li ne ratio, but likely requires a low in the active 
region given the low observed ratios of to intensity in IWeifi et al.l (j2QQ5l ). 

A more realistic, but more complex, prescription adopts ^ 2.5 Mq pc“^ (K km s“^)“^ in the extended disk 

and the faint part of the outflow and takes Gqq^ I Mq pc“^ (K km s“^)“^ in starburst region and the bright part 
of the outflow, which is still very near the disk. A convenient cut for these purposes is Iqo ~ 20 K km s~^ in the 
outflow region: 


r2.5 Mq pc ^ (K km s extended disk, elsewhere if Ico < 5 K km s 

\l.O Mq pc“^ (K km starburst, outflow where Ico > 5 K km s“^. 


(B8) 
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This approach posits that CO is a more effective emitter in the bright starburst region than the cooler, more quiescent 
extended disk. This agrees with the sub-mm and dust fitting results well. However, with coarse resolution and only a 
single line for the very extended part of the outfiow, we lack concrete knowledge of what the CO emitting structures in 
the wind look like in detail. Bear in mind that this two-regime prescription is specific to our data because we expect 
the intensity to be a strong funct ion of angular resolution. 

Figure [24] shows that Equation IBS! yields results that are consistent with a single, linear relation between dust and 
gas over a large area in M82 a nd that this approach re presents a substantial improvement over using only the typical 
value for disk galaxies found bv iSandstrom et all (j2Q13f ). The figures show that the histogram of dust to gas ratios over 
and individual line of sight and the binned scaling relation between dust and gas surface density both become more 
consistent with a single dust-to-gas ratio using this approach. Note that the bright starburst itself does not appear to 
follow this prescription. Here, the conversion factor may become even more complicated as it convolves many distinct 
sets of physical conditions a nd o pacity becomes important; however, the burst itself is not the topic of this paper. For 
the rest of galaxy. Equation IB8I results in single linear scaling between dust and gas. 

TESTS RELATED TO THE DUST TREATMENT 

Figure [25] illustrates fo ur tests related to o ur treatment of IR emission around M82. The top two panels show the 
effect of our choice of the iDraine T3 (l2nnl model by comparing these to the result of a simple modified blackbody 
fit to the IR SED between 70 and 350/im. The top left panel shows our modeled Edust as a function of the optical 
depth at 250/im from the blackbody fit with each pixel as a single gray point. The two track each other closely, 
with about a factor of two scatter and a similar systematic trend across the range of dust surface densities in the 
galaxy. The normalization needed to equate the two would correspond to a mass absorp tion coefficient at 25Q/im of 
K ^ 1-2 cm^ g~ ^ . Th is is moderately low compared, e.g., to the theoretical value (e.g., IDraine fc Lmnni because 
the IDraine fc M il2nnl models model a population of dust and so tend to yield somewhat higher dust mass surface 
densities than single population fits. The treatment of heating also differs and matters to the derived surface density 
(the right panel). Given that most of the results in this paper are comparative and do not depend on the normalization 
of the dust models Figure [25] shows that our adopted model has only a modest impact. Th e right panel compp es the 
derived temperature from the modified blackbody fit and the mean radiation field from the IDraine fc Lil (j2QQ7f ) model 
fit. For comparison the line expec ted if T = 20 K cor responds to U = 1 and jS = 1 appears as a black line. The range 
of radiation fields modeled by the IDraine fc Lil (j2QQ7[ ) model lead to a lower mean U at low dust temperatures, which 
is a goal of the model and a more reali stic treatment of the SED for low resolution observations of other galaxies (see 
iDale fc Heloull2QQ2l: IDraine et al.ll2QQ7h . At high temperature, high radiation field regions, the two treatments tend to 
converge. The two quantities do track one another, though not linearly. Given that use of {U) is mostly qualitative, 
the comparison shows that our results for {U) could mostly be rec ast in terms of dust temperature if one prefers to 
avoid the assumptions about unresolved populations underlying the IDraine fc Lil ()2QQ7[ ) model. 

The bottom left panel illustrates the effect of subtracting a median background from the SPIRE maps when con¬ 
structing the large-scale profiles. We plot the profiles with and without background subtraction, illustrating that the 
wide maps around M82 do converge to a non-zero level at large radii. As seen in the main text, this emission does 
not appear to be highly structured or to resemble emission at other wavelengths and the Herschel data processing is 
not intended to provide a reliable zero-point intensity. The figure illustrates a plausible range of background levels via 
the shaded regions and then shows the effect of varying the subtraction within this range on the subtracted profiles. 
Within the ~ 8 kpc from the galaxy center this treatment appears fairly robust. The nature of faint, extended emission 
and the presence of any faint, smooth emission on very large scales around M82 will have to be addressed with facilities 
or data processing optimized for that purpose. 

Finally, the bo t tom r ight panel of the plot the expected ratio of collision to radiative heating for a simple treatment 
following iDraind (j2Qllf ) . They provide an estimate of this ratio for atomic hydrogen, their Equation 24.7. We show 
the results of this calculation for a range of densities and temperatures, a solar neighborhood radiation field {U = 1), 
which is already very low given our results on M82. We adopt a grain size 0.1/im and an equal mix of graphite and 
silicate grains following their treatment; larger grains would only increase radiative heating. The plot shows that for 
any plausible set of conditions in the bulk of the molecular or atomic gas, radiative heating is expected to dominate 
over collisional heating in M82. Shocks have been shown to be important in a variety of ways to both the dust and 
the gas around M82 (e.g., iKaneda et al.ll2QT(l iBeirao et al.ll2Ql^ but we do not expect that they push the bulk of the 
gas mass into the extreme regime needed for collisional heating to dominate radiative heating in the presence of the 
strong radiation fields around M82. In individual shocked regions or the dense, well-shielded parts of the starburst, 
the picture may be different. For the large-scale view we study, where neutral gas and dust show good correspondence, 
we know that the bulk H 2 temperature is likely a few times 10 K and constrained by IR spectroscopy to be not 
more than a few times 100 K . Meanwhile, the hot gas at extreme temperatures is expected to have low densities 
(|Strickland fc Heckmanll2009f ). even though its temperature is extreme. More detailed study is needed to understand 
dust heating in the interface of the hot wind material and the cold material, but this first order calculation supports 
the interpretation of dust temperature as a tracer of the radiation field around M82. 

THE EEEECTS OE THE IRAM 30-M ERROR BEAM 

M82 represents a prime example of the general astronomical problem of faint emission in the presence of a bright 
source. Although the problem is far less several in mm-wave line emission than at other wavelengths, stray light pickup 
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Fig . 25.— Calculations related to our dust treatment. Top row: comparisons of modeling using the IDraine &: Lil (|2003)JDrainee^^lJ 
(|2007h models to fitting a single modified blackbody to the infrared SED. (top left) The dust surface density from the IDraine fc Lil (|2007l) 
modeling as a function of modified blackbody optical depth, a tracer of dust column. The two track one another well mo dulo a normal i zation 
factor (the dashed lines are separated by a factor of two with slope unity), (top right) Mean radiation field from the IDraine &: 
modeling as a function of dust temperature from a blackbody fit with a simple expectation as a solid black line. The two converge in 
the hi gh-temperature, h^ h radiation field regime. They track one another almost monotonically at lower temperatures but differ because 
of the IDraine fc Lil (120071 ) treatment of an unresolved population. The comparisons show that the results of the paper would not change 
substantially if we shifted our dust treatment, (bottom left) The need for and effect of the subtraction of a median background from the 
SPIRE maps. We show averaged annular profiles around M82 at 250 and 350/rm before and after subtraction of a median background. We 
also indicate a plausible range of zero-point values (light rectangular shaded regions) and the effect of varying the background within this 
range on the subtracted profiles (darker shaded region). The effect is minimal within the region that we study in t he paper (out to the 
black line), (bottom right) Simple calculation of the balance of collisional and radiative heating following [Drm^ (| 201 H i for a mean radiation 
field U = 1. Radiative heating dominates collisional heating for the plausible range of bulk gas temperatures and densities, supporting 
the interpretation of the dust temperature in terms of illuminating radiation field, though more detailed modeling of the interface of cold, 
dense and hot, rare material is needed. 


can still be a problem. The main beam efficiency of the IRAM 30-m, 77nib, is measured to be ~ 58% at 230 GHz. 
Much of the rest of the response to emission on the sky is in an “error beam” that has been characterized as a series 
of Gaussians much larger in size than the main beam. As a result, our map could include low-level pickup of extended 
structures. In fact, the large size of the error beam and the fact that it operates “by channel” (i.e., at a specific 
frequency) combine to make this a small effect in M82 and indeed throughout the HERAGLES survey. This appendix 
describes the quantitative test that we use to arrive at this conclusion. 

Gonsider a data cube that has been reduced and converted to an image on the sky under the assumption the error 
beam is free of astronomical signal, so that the cube in units has been converted to T^b units by multiplying by a 
factor ? 7 “^. This cube represents an upper limit to the true sky emission because some emission observed with the error 
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beam may be included in the cube. Under the assumption that the error beam effect is a perturbation on the overall 
emission, we can estimate the magnitude of such contamination by simply convolving the observed cube with the error 
beam plane by plane and then scaling the result by (we scale to account for the fact that this multiplication by 
^mb already applied to the contamination in the original cube). This now represents an estimate of the error beam 
contamination in the Tmb cube that can be subtracted from the original cube. Often this calculation will be enough to 
demonstrate the minimal impact of the error beam. In the case where the contamination is significant, the calculation 
may need to be refined because this first estimate of contamination will be an overestimate (recall that the original 
cube itself was an overestimate because it includes the contamination). To refine the estimate, one subtracts the first 
estimate of error beam contamination from the original cube and then convolves the new, contamination-corrected 
cube with the error beam. This is a refined estimate of contamination by the error beam that can be subtracted from 
the original cube to yield an improved estimate of the true intensity distribution. As long as the error beam represents 
a modest correction to the overall emission, iteration should yield a corrected estimate of the true emission. 

Procedurally, the correction scheme is: 

1. Generate a cube in Tmb units assuming no contamination. 

2 . Convolve the cube with the known error beam plane by plane. Scale by 77 “^. This represents an estimate of 
error beam contamination. 

3. Subtract the estimated contamination from the original cube. This is an estimate of the true, error-beam 
corrected sky distribution in Tmb units. 

4. Return to step #2 and use the corrected cube to generate the contamination estimate. 

5. Repeat steps #2-4 until the change between successive iterations is small. 

In practice, even this iterative approach often yields an overestimate of the effect of the error beam on single dish 
maps of other galaxies because the procedure of baseline fitting will tend to suppress stray light pickup in the original 
cube. The error beam is often comparable to the size of the galaxy, so that its effect is to place a heavily reduced 
version of the spectrum of the galaxy into individual spectra. For a highly resolved galaxy, this spectrum will often 
be quite extended compared to the line profile at a given location within th e galaxy. Most baseline fitting procedures 
(for example the local Hi and CO windowing function in lLerov et al .11200^ will remove such l ow, broad emis s ion. A 
corollary of this, the absence of broad line wings in stacked HERACLES spectra, is discussed in iSchrnba et al.l (j2Qllf ). 

In an IRAM memo, Kramer, Pehalver, and Greve (2013) provide the latest characterization of the IRAM 30-m 
beam, the forward efficiency, and forward spillover efficiency term that is the difference between the main beam and 
forward efficiencies. Adopting their prescription for the error beam, we iteratively model the error beam for M82 (and 
all of the other HERACLES) targets following the procedure above. We then re-fit spectral baselines in the cubes 
following the same robust polynomial fitting procedure used to produce the original cubes. We iterate a total of five 
times for M82, though the last three iterations have only small effect. The result of this modeling for M82, which 
should be among the worst cases for the HERACLES galaxies, is a decrease of ^ 10% in the overall flux of the cube. 
Locally, the effect can be more extreme, with ^ 20-30% differences between the corrected and uncorrected maps in 
the area immediately around the bright starburst disk. 







